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ABSTRACT 

We consider the possibihty that at least some GRB explosions take place inside 
pulsar wind bubbles (PWBs), in the context of the supranova model, where initially a 
supernova explosion takes place, leaving behind a supra-massive neutron star (SMNS), 
which loses its rotational energy on a time scale of months to tens of years and collapses 
to a black hole, triggering a GRB explosion. The most natural mechanism by which 
the SMNS can lose its rotational energy is through a strong pulsar type wind, between 
the supernova and the GRB events, which is expected to create a PWB. We analyze 
in some detail the observational implications of such a plerionic environment on the 
afterglow and prompt GRB emissions, as well as the prospect for direct detection of 
the plerion emission. We find that for a simple spherical model, GRBs with iron lines 
detected in their X-ray afterglow should not have a detectable radio afterglow, and 
should have small jet break times and non-relativistic transition times, in disagreement 
with observations for some of the GRBs with X-ray lines. These discrepancies with 
the observations may be reconciled by resorting to a non-spherical geometry, where the 
PWB is elongated along the polar axis. We find that the emission from the PWB should 
persist well into the afterglow, and the lack of detection of such a component provides 
interesting constraints on the model parameters. Finally, we predict that the inverse 
Compton upscattering of the PWB photons by the relativistic electrons of the afterglow 
(external Compton, EC) should lead to high energy emission during the early afterglow 
that might explain the GeV photons detected by EGRET for a few GRBs, and should 
be detectable by future missions such as GLAST. 

Subject headings: gamma rays: bursts — pulsars: general — supernova remnants — stars: 
winds, outflows — shock waves — radiation mechanisms: nonthermal 



1. Introduction 

The leading models for Gamma-Ray Bursts (GRBs) involve a relativistic wind emanating 
from a compact central source. The prompt gamma-ray emission is usually attributed to energy 
dissipation within the outflow itself, due to internal shocks within the flow that arise from variability 
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in its Lorentz factor, while the afterglow emission arises from an external shock that is driven into 
the ambient medium, as it decelerates the ejected matter (Rees & Meszaros 1994; Sari & Piran 
1997). In this so called 'internal-external' shock model, the duration of the prompt GRB is directly 
related to the time during which the central source is active. The most popular emission mechanism 
is synchrotron radiation from relativistic electrons accelerated in the shocks, that radiate in the 
strong magnetic fields (close to equipartition values) within the shocked plasma. An additional 
radiation mechanism that may also play some role is synchrotron self-Compton (SSC), which is the 
upscattering of the synchrotron photons by the relativistic electrons, to much higher energies. 

Progenitor models of GRBs are divided into two main categories. The first category involves 
the merger of a binary system of compact objects, such as a double neutron star (NS-NS, Eichler 
et al. 1989), a neutron star and a black hole (NS-BH, Narayan, Pacyhski & Piran 1992) or a 
black hole and a Helium star or a white dwarf (BH-He, BH-WD, Fryer &; Woosley 1998; Fryer, 
Woosley & Hartmann 1999). The second category involves the death of a massive star. It includes 
the failed supernova (Woosley 1993) or hypernova (Pacyhski 1998) models, where a black hole is 
created promptly, and a large accretion rate from a surrounding accretion disk (or torus) feeds a 
strong relativistic jet in the polar regions. This type of model is known as the collapsar model. An 
alternative model within this second category is the supranova model (Victri & Stella 1998), where 
a massive star explodes in a supernova and leaves behind a supra-massive neutron star (SMNS) 
which on a time scale of a few years loses its rotational energy and collapses to a black hole, 
triggering the GRB event. Long GRBs (with a duration > 2 s) are usually attributed to the second 
category of progenitors, while short GRBs are attributed to the first category. In all the different 
scenarios mentioned above, the final stage of the process consists of a newly formed black hole with 
a large accretion rate from a surrounding torus, and involve a similar energy budget (< 10^^ ergs). 

In this work we perform a detailed analysis of the supranova model, focusing on its possible 
observational signatures. This aims towards establishing tools that would enable us to distinguish 
between the supranova model and other progenitor models through observations, and to constrain 
the model parameter using current observations. The original motivation for the supranova model 
was to provide a relatively baryon clean environment for the GRB jet. As it turned out, it also 
seemed to naturally accommodate the later detection of iron lines in several X-ray afterglows 
(Lazzati, Campana, & Ghisellini 1999; Piro et al. 2000; Vietri et al. 2001). 

It was later suggested that the most natural mechanism by which the SMNS can lose its 
rotational energy is through a strong pulsar type wind, between the supernova and the GRB events, 
which typically creates a pulsar wind bubble (PWB), also referred to as a plerion (Konigl &: Granot 
2002, KG hereafter; Inoue, Guetta &: Pacini 2002, IGP hereafter). KG suggested that the shocked 
pulsar wind into which the afterglow shock propagates in this picture may naturally account for 
the large inferred values of eg ~ 0.1 and ~ 10~^ — 0.1 (the fractions of the internal energy 
in the electrons and in the magnetic field, respectively) that are inferred from fits to afterglow 
observations (Wijers & Galama 1999; Granot, Piran & Sari 1999; Chevalier & Li 2000; Panaitescu 
&L Kumar 2002). This is attributed to the fact that pulsar winds are believed to largely consist 



-3- 



of electron-positron pairs, and have magnetization parameters in the right range. This relaxes the 
need of generating strong magnetic fields in the shock itself, as is required in other models, where 
the magnetic field in the external medium (assumed to be either the ISM or a stellar wind of a 
massive star progenitor) is typically too small to account for the values of that are inferred from 
observations. Another attractive feature of this model, pointed out by IGP is the possible high 
energy emission, in the GeV-TeV range, that may result from the upscattering of photons from 
the plerion by the rclativistic electrons in the afterglow shock (external Compton, EC hereafter), 
and may be detected by GLAST. They have shown that the EC emission can provide a viable 
explanation for the extended GeV emission seen by EGRET in GRB 940217 (Hurley et al. 1994). 

We use a simple spherical model for the PWB. We find that a spherical model cannot accom- 
modate the typical afterglow emission together with the iron line features observed in the X-ray 
afterglow of some GRBs. However, it was mentioned early on that in order to have a long lived 
afterglow emission together with the iron line features, a deviation from spherical symmetry is 
needed, where the line of sight is relatively devoid of the material producing the iron lines (Lazzati 
et al. 1999; Vietri et al. 2001). This is required in order to avoid a direct collision of the afterglow 
shock with the line producing material on an observed time of the order of a day or so. It was later 
pointed out that a PWB is expected to exist inside the SNR shell, which decelerates the afterglow 
shock at a smaller radius, so that in order for the afterglow to remain rclativistic up to a month or 
more, and produce the iron lines, we need the PWB to be elongated along its rotational axis (KG). 
In this paper we strengthen this conclusion, and show that in order to produce iron lines with a 
spherical PWB, its radius must be sufficiently small, resulting in a large density inside the PWB 
and a high self absorption frequency implying no radio afterglow, in contrast with observations. 
We leave the detailed treatment of an elongated PWB to a future work, while in the present work 
we briefly comment about the expected effects of an elongated geometry compared to a spherical 
one. 

In this work we extend the analysis of KG and IGP, and perform detailed calculations of the 
radiation from the PWB, the prompt GRB and from the afterglow that occur inside the PWB. 
We now give a short overview of the structure of the paper, where in each section we stress the 
original features, new results and the observational constraints on the model. In §2 wc present 
our "PWB" model, introduce the relevant parameterization and model the acceleration of the 
supernova remnant (SNR) shell by the shocked pulsar wind. We use a simple spherical geometry 
and the pulsar wind is assumed to consist of proton and components with roughly equal energies, 
as well as a magnetic field. The conditions under which the iron line features that were observed 
in several X-ray afterglows may be reproduced within the PWB model, are investigated in §3. We 
find that this requires a time delay of < 1 yr between the supernova and the GRB events. In §4 we 
perform a detailed study of the plerion emission, including the synchrotron and SSC components, 
and provide an elaborate description of the relevant Klein-Nishina effect. Wc also discuss the upper 
cutoff that is imposed on high energy photons due to pair production with the radiation field of the 
PWB, go over the prospect for direct detection of the plerion emission, and derive observational 
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constraints on the parameters of our model. The effects of the PWB environment on the prompt 
GRB emission are analyzed in §5, and we find that the EC from the prompt GRB should typically be 
very small, but might be detectable for extreme parameters. In §6 wc discuss the implications of a 
plerionic environment on the afterglow emission, and introduce the appropriate parameterization. 
The radial density profile of the PWB is approximated as a power law in radius (KG), oc r~*^, 
where k typically ranges between (similar to an ISM) and 1 (intermediate between an ISM and 
a stellar wind). The synchrotron, SSC and EC components are calculated and wc provide detailed 
expressions for the break frequencies and flux normalization, for /c = 0, 1. We also calculate the 
high energy emission that is predicted in this model. The results are discussed in §7 and in §8 we 
give our conclusions. 



2. The Pulsar Wind Bubble 

Within the framework of the supranova model, a SMNS (also simply referred to as a pulsar) 
is formed in a supranova explosion, and then loses a large part of its rotational energy before 
collapsing to a black hole and triggering the GRB. The most plausible mechanism for this energy 
loss is through a pulsar type wind (KG; IGP). A pulsar wind bubble (PWB) is formed when the 
relativistic wind (consisting of rclativistic particles and magnetic fields) that emanates from a pulsar 
is abruptly decelerated (typically, to a Newtonian velocity) in a strong rclativistic shock, due to 
interaction with the ambient medium. When a bubble of this type expands inside a SNR, it gives 
rise to a plerionic SNR, of which the Crab and Vela remnants are prime examples. Motivated by 
previous works (Rees & Gunn 1974; Kennel k. Coroniti 1984; Emmering k, Chevalier 1987; KG) we 
consider in detail a spherical model where the shocked pulsar wind remains largely confined within 
the SNR. In §7 briefly discuss the possible consequences of some more complicated geometries. 

The wind luminosity may be estimated by the magnetic dipole formula (Pacini 1967), 

^- = -6^ = '-°"^° [to^) (isk^J [i6^) ' 

where is the polar surface magnetic field, i?* is the circumferential radius (neglecting the dis- 
tinction between its equatorial and polar values in this approximation), and ft^ is the (uniform) 
angular velocity (whose maximum value is ~ 2 x 10^ s~^; e.g., Haensel, Lasota, &; Zdunik 1999). 
The spin-down time of a rapidly rotating SMNS can be estimated as 

*^'^= l7"H^](,5:^J [t5^) [w^y [w^J ^ ^ 

(see Vietri &: Stella 1998), where -Erot = aGM^Q,^/2c is the portion of the rotational energy of 
an SMNS of mass and angular velocity 0^, that needs to be lost before it becomes unstable 
to collapse.^ The spindown timescale, that sets the time delay between the supernova and GRB 



^The total rotational energy of the SMNS is given by jGM»fl*/2c, where the parameter j measures the stellar 
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events, depends on the physical parameters of the SMNS. Of these parameters, the least constrained 
is the magnetic field, which is typically expected to be in the range of ~ 10"*^^ — 10"*^'^ G, and may 
cause a variation of > 2 orders of magnitude in tsd- There is also a strong dependence on the radius 
i?*, which depends on its mass and the (uncertain) equation of state, that may account for a 
change of up to ~ 1 order of magnitude in the scaling of tgd- For example, for i?* = 10 Km, with 
the values of the other parameters as given in Eq. (2), we have tgd ~ 60 yr. We conclude that the 
expected range of tgd is from a few weeks to several years. 

During tgd the luminosity of the wind is roughly constant and the wind should energize the 
PWB depositing an energy of the order of E^ot- The luminosity of the pulsar wind is divided 
between its different components: fractions ^pj and a^, in pairs, protons and Poynting flux 
(magnetic field), respectively. The inferred values of = (^w/i^ + (^w) for PWBs, such as Vela 
or the Crab, are typically ~ 10~^ (Arons 2002), though there are also estimates as high as ~ 1 
(Helfand, Gotthelf Sz Halpern 2001). We shall adopt a fiducial value of = 10~^, which implies 
+ = 1 ~ — 1- Gallant & Arons (1994) inferred ^p/^e ~ 2 for the the Crab, and we adopt 
this estimate for our fiducial values, and use = 1/3 and = 2/3. The inferred values of the 
Lorentz factor of the pulsar wind are ~ 10^ — 10^. 

The SN ejecta is accelerated by the force exerted due to the pressure of the expanding PWB, 
Poutj at its outer boundary, i.e. at the radius of the SNR, i?sNR- The pressure is expected to drop 
by a factor of order unity between the radius of the wind termination shock, Rg, and i?sNR (KG), 
so that Pout ^ VpPcLVj 

where pav = V^w/'^y = rjEyj/AT^B? is the average pressure inside the PWB, 
r] is the fraction of the energy E^q^ that remains in the PWB, V = (47r/3)i?3 is the volume of the 
PWB (which can be different for an elongated PWB) and E^ Erot{t/tsd) is the energy emitted 
in the wind up to the time t. As we show below, the PWB is typically in fast cooling, and the 
electrons lose all of their internal energy to radiation, implying ^7 ~ + ~ ^p- The equation of 
motion, Msnr-R = 47ri2^Pout! may be written as 

RR=^f^t. (3) 



In the case of a non-spherical PWB, MgNR should be replaced by the isotropic equivalent mass, 
-^iso(^) = 47rdMsNR/(^^^, and R becomes R{9), where 6 is the angle from the polar axis. However, 
the pressure can be taken as independent of the angle 6, since the shocked pulsar wind is highly 
sub-sonic. If MisoiO) is smaller near the poles {6 ~ 0) and larger near the equator (6 ~ vr/2), then 
even if the SNR shell is initially spherical, the acceleration would be larger near the poles, resulting 
in a much larger polar radius, Rp, compared to the equatorial radius, R^q. Rp ^ Req. This is a 
natural mechanism that can lead to an elongated geometry for the PWB. 

Returning to the spherical case, the acceleration becomes significant at the time iacc (and 



angular momentum in units of GM^/c and has values in the range 0.57 — 0.78 for realistic equations of state (e.g., 
Cook, Shapiro & Teukolsky 1994; Salgado et al. 1994). 
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radius i?acc ~ 'I'o^acc) when Eyj first exceed the initial kinetic energy of the SNR, Eq = Ms-nrVq/2 
(~ 10^^ ergs), where tacc = isd-^oZ-Srot • We therefore have 

tsd E^ot fvbV f^^V^'^wO, (4) 



''acc E!o V ^0 / \ -Race 

where Vb = ^^SNR(isd) smd Rb = -RsNR(isd)- The dynamics are given by 

Vot t < tacc 



i-acc 

Vo t < tacc 



(5) 



(6) 



''acc 



These scahngs cigree with the results of Reynolds Chevalier (1984). At t ^acc 

we have 

Esnr/Ew = 3ripri/2, and conservation of energy imphes that Eq + E^ « E^ = rjEyj + i?sNR = 
(3r/p/2 + l)'qEw or (3??p/2 + 1)77 = 1. Our fiducial value of = 2/3 implies r^p = 1/3, which 

is reasonable. We also obtain that -EsNR(^sd) = (3??/2)£^rot ~ -E'rot- For a typical ejected mass, 
-^SNR ~ IOM0, this would imply ~ 0.1c. Finally we have 

Rb = \vbt,d = 6.3 X 10iV6,-iisd,o (7) 

where we set V},/c = Pj, = ^-^Pb-i and tgd = isd,o y. To the extent that Vh oc (-Ej-ot/MsNR)^''^ has 
nearly the same value in all sources, the magnitude of Rf, is determined by that of tgd- In a similar 

vien, if the energy lost during the SMNS lifetime, ^^j-ot = 10^^ £'53 ergs, is approximately constant 
from source to source (£'53 ~ 1), then tsd can also be used to parameterize the SMNS wind power: 
L„ = £rot/tsd = 3.2 X 1045£53/ty, ergs s-\ 

The acceleration of the SNR shell by the lower-density bubble gas would subject it to a 
Rayleigh- Taylor (RT) instability, which could lead to clumping (Jun 1998). The growth timescale 
of the RT instability on a spatial scale R is tRT ~ {R/RY^^. The important quantity to estimate in 
order to see if the SNR shell can be clumped is the ratio between tux and the dynamical timescale 

tdyn = (R/R) 

tRT _ R _ / 2£'sNR /„^ 

(it;i?)V2 Y ripTjEyj ' 

where we have used Eq. (3). This implies tRx/iexp ~ during the acceleration (tacc < t < tgd)- 
This could produce only moderately strong fragmentation over the dynamical time of the system. 
However, as the acceleration occurs over ~ 3 orders of magnitude in radius [see Eq. (4)], the radius 
doubles itself ~ 10 times during the course of the SNR acceleration, so that despite the fact that 
tRT I texp is of order unity, it is feasible that considerable dumpiness may still be caused due to the 
RT instability. An even stronger fragmentation may occur if the RT instabilities grow on a length 
scale X ~ aR smaller than R (i.e. a < 1), where in this case tRx/iexp ~ \/3a. 
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In order to calculate the emission from the plerion, we use the average quantities of the shocked 
pulsar wind within the PWB, and neglect their variation with radius. The latter is expected to be 
more important in the afterglow emission, and is therefore taken into account in §6 that discusses 
the afterglow emission. The postshock energy density is given by 

where V = (47r/3)i?| is the total volume within the PWB, and is approximately equal to the volume 
occupied by the shocked wind. For an elongated PWB the expression for the volume, V, will be 
different, but it can directly be plugged into these equations, in place of the spherical expression. 
The injection rates of electron-positron pairs and of protons at the source are given by 

iv,,, = . (10) 

Hence, the total number of particles within ii^ at time t is Ne^p{t) = N^^pt, and the number density 
at is 

^ -^e,p(^sd) Ce,p 6 2^e,p-^rot (11) 

Fractions e(,B, e{,e, e^bp of the post-shock energy density go to the magnetic field, the electrons and 
the protons, respectively. We expect these fractions to be similar to those in the pulsar wind 
(cfeS ~ ~ "^loi ^be ^ ie) and usc the corresponding fiducial values. Subscripts containing the 
letter 'b' denote quantities related to the PWB. The electrons will lose energy through synchrotron 
emission and inverse-Compton (IC) scattering. We study the characteristic features of the plerion 
emission and investigate the conditions required for the production of the observed iron lines and 
whether the plerion can be detected by the present instruments. Moreover, an important implication 
of the plerion emission is that the GRB should explode inside a radiation rich environment (i.e., 
the luminous radiation field of the PWB). The external photons arc highly Dopplcr-boostcd in 
the rest frame of the shocked fluid, for both internal and external shocks (that are responsible for 
the prompt GRB and afterglow emission, respectively), and can act as efficient seed photons for 
IC scattering (external Comptonization, EC). We study the observational consequences of the EC 
process, both for the prompt GRB emission and for the afterglow. 

Since we use a large number of parameters in the paper, and in order to make it easier to follow 
all the different parameters, we include a table (Table 1) with the most often used parameters, where 
we mention the meaning of each parameter and the fiducial value that we use for that parameter. 



3. X-ray Lines in the Afterglow 



One of the main motivations for the supranova model is that it can naturally explain the 
detections of iron lines in the X-ray afterglows of several GRBs, both in emission (GRB 970508, 
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Table 1. The parameters most often used in the paper, their meaning and fiducial values. 



Parameter 


meaning 


fiducial value 




time delay between SN and GRB 


1 or 10^-^ yr 




Lorentz factor of pulsar wind 


104.5 


Erot 


rotational energy lost by SMNS 


10^3 erg 


MsNR 


the mass of the SNR shell 


10 Mq 




fraction of E'rot that remains in the PWB 


2/3 




fraction of wind energy in magnetic field 


10-3 




traction oi wmd energy in e pairs 


1 /o 

1/3 




fraction of wind energy in protons 


2/3 




velocity of SNR shell at tsd (in units of c) 


0.1 


^bB 


fraction of PWB energy in magnetic field 


10-3 


^be 


fraction of PWB energy in pairs 


1/3 


a 


fraction of e"^ energy that is radiated 


1 


S 


power-law index of distribution in PWB 


2.2 


P 


electron power-law index in the GRB 


2.5 


r 


bulk Lorentz factor of GRB ejecta 


102.5 




kinetic luminosity of GRB outflow 


1052 erg/s 


ty 


variability time of GRB central engine 


10 ms 


k 


n, e oc r~'^ in the PWB 


or 1 




fraction of AG energy in magnetic field 


10-3 




fraction of AG energy in electrons 


0.1 


Eiso 


isotropic equivalent energy in AG shock 


1053 erg 


t 


observed time since GRB 


1 day 



-9- 



Piro et al. 1998; GRB 970828, Yoshida et al. 2001; GRB 000214, Antonelli et al. 2000; GRB 

991216, Piro et al. 2000) and in absorption (GRB 990705; Amati et al. 2000). The statistical 
significance of these detections is at the ^ 3a level, with the exception of GRB 991216 where a k-a 
emission line was detected with a significance of ~ 4c7. Emission lines of lighter elements (Mg, Si, S, 
Ar, Ca) have been reported in the X-ray afterglow of GRB 011211, at the level of 3cr (Reeves et al. 
2002). This latter detection has been disputed by other authors (Borozdin &; Trudolyubov 2002; 
Rutledge &; Sako 2002), and may be said to be controversial. These line features may naturally arise 
in the context of the supranova model, where a SNR shell is located at a distance of R> 10^^ cm 
from the location of the GRB explosion (Lazzati, Campana, & Ghisellini 1999; Piro et al. 2000; 
Lazzati et al. 2001; Vietri et al. 2001; Bottcher, Fryer, &; Dermer 2001). In this section we explore 
the condition under which such features may occur within our model, and obtain the relevant 
constraints on the model parameters. In the following sections we investigate the implications of 
these constraints on the other observational signatures of the model: the plerion, prompt GRB and 
afterglow emissions. 

We derive constraints on our model parameters using the observational data for GRB 991216, 
as an example of an afterglow for which iron lines were detected, since this is the most statistically 
significant detection to date. Similar constraints may be obtained for other afterglows with X-ray 
features, using similar arguments. The X-ray afterglow of GRB 991216 was observed by Chandra 
from 37 hr to 40 hr after the GRB, and shows an emission line at 3.49ib0.06 KeV, with a significance 
of ~ 4cr. The line flux was Fpe ~ 1.6x10"^^ erg cm~^ s~^, which for a redshift z = 1.02 of this burst 
(Vreeswijk ct al. 2000) implies an emission rate iVpc 4 x 10^^ of line photons, a luminosity of 
Lpe ~ 4 X lO'^^erg and a total energy of Epc 3 x 10'^^ ergs assuming the line emission lasted 
for tpe ~ 40/(1 -I- z) hr in the cosmological rest frame of the GRB (Vietri et al. 2001). 

In the simplest version of our model, we assume a spherical geometry, and identify the line 
emitting material with the SNR shell, that is located at a radius Rf,- We use the above observations 
to derive constraints on R^, or equivalently, on the time delay between the SN and the GRB events, 
tsd- The value of Rb may be constrained by the requirement that the geometrical time delay in the 
arrival of the photons to the observer, ~ Rij9^^^/2, should not exceed the total duration of the iron 
line emission, tpe, 

«'S1.7xlO-(^)(^)"cn, . ,,<27(|i^)(^)"ft- ^, (12) 

where we have identified the opening angle to which the ionizing radiation extends, ^j-ad) with the 
jet opening angle, Oj 0.05 (Frail et al. 2001; Panaitescu & Kumar 2002). 

Another constraint may a^risc from tli6 rcQiiircmcnt that A^Pe — -^Fe/^rcc ^ 

4 X 10^^ s \ where 

the recombination time is given by Uec ~ 4 x 10^ Z-'^T^-^n-^ = 2.8 x lO^^TgO-^n-^ s, Tg = Tg/lO^ K, 
and we have assumed an electric charge of -Z^ = 24 for the iron ions (Lazzati et al. 2001; KG). In 
order to parameterize the electron number density, rig, we need to relate between the width of the 
SNR shell, ARf,, and its radius, Rj,. If the SNR shell is efl&ciently fragmented during its acceleration 
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phase, due to the RT instability, then one might expect dense clumps of size Id <C Ai?f,, spread over 
a radial interval of ARf,, that cover a fraction of order unity of the total solid angle. This amounts 
to an effective width for the SNR shell of Ai?e// = Id ^ ^Rb, implying rig = MsNR/47ri?^Ai?e//"^p 
and 

R,<9x W''M's{,l^,Mi^_,Cl^%'/' cm , t^^ < lAM'si,n,,Mi^_,Cl^%'^' P^,'^, yr , (13) 

where ^ = AReff/Rb = 10~^^_3 and Mpe = O-lMpe-iM© is the mass of the iron in the SNR shell. 
The SNR shell is compressed during its acceleration, and may attain AR^ < OARj,, so that ^ may 
be as low as ^ 10~^. 

A final constraint may be derived by considerations related to the total energy budget. The 
total energy in the line is E-pc = eEj ~ 3 x 10^^ ergs where the efficiency e is the product of the 
ratio of energies in the ionizing X-ray continuum and the prompt gamma-ray emission, and the 
energy fraction of the X-ray continuum that goes into the line emission, and is expected to be 
< 0.01 (Ghisellini et al. 2002). This implies > 3 x 10^^ ergs which is somewhat in excess of 
the value E^ ^ 7 x 10^*^ ergs found by Frail et al. (2001). If the optical depth of the iron atoms 
is Tpe < 1, then the efficiency is further reduced by a factor of rpe = MpQaFe/ATrR^SSnip, where 
(Tpe ~ 2.0 X 10~^° cm^ (Krolik & Kallman 1987). Therefore, we must have rpe > 1, i-e. 

i?6 < 6 X lO^^M^/li cm , tsd < 0.9M^/li/3,;^ yr , (14) 

We conclude that X-ray line features, such as the ones observed in the afterglow of several 
GRBs, may be accommodate in a spherical model only for tgd ^ 1 yr, Rb ^ 10^'' cm. This 
constraint can be relaxed for an elongated PWB. In this case the iron lines may be produced by 
the material near the equator, which is at a much smaller radius than the polar radius, enabling 
the afterglow shock, that propagates along the poles, to reach a considerably larger radius. 



4. The Plerion Emission 

In this section we evaluate the luminosity and the spectrum emitted by the plerion. We use 
the average values of the quantities within the PWB, that were derived in §2. As we show below, 
the electrons are in the fast cooling regime for relevant values of tsd, and therefore most of the 
emission takes place within a small radial interval just behind the wind termination shock, and the 
various quantities should not vary significantly within this region, and should not be very different 
from their average values within the PWB. 

The magnetic field inside the plerion is 

Bb = VSTTCbBe ^ 1.3 V2/lHB,-3^lil53l3b,~lhdfi ^ ' (1^) 

where ebB,-3 = ebs/lO^ and 772/3 = 77/(2/3). Relativistic electrons/positrons (hereafter simply 
electrons) are injected into the plerion at the rate Ng (see Eq.lO)with a power law distribution 
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iV(76) = dn/djf, oc 7j * in the Lorentz factor (LF) range 7;,^ ^ 7& ^ 76M- The minimum electron 
LF is given by 

(s — 2\ eft e 
7^ ) '7^^ ~ ^'^ ^ 10^%/3e6e,l/3^~V3T'"'.4-5 > (16) 

where 7^^,4.5 = 7«)/10^-^, (-he.ilz = e6e/(l/3), ^e,i/3 = Ce/(l/3) and we use s = 2.2 to obtain the 
numerical values for the rest of the paper. The electrons radiatively cool by the combination of 

the synchrotron and synchrotron-self-Compton (SSC) process, the timescales of which arc tgyn ~ 
Q-Kruec/ gtB'^^i, and tsc = tsyn/Yb and the combined cooling time being tc = (l/tgyn + ^/tsc)~^ = 
*syn/(l + Yb), where 

a^be/^bB aebe/^bB ^ 1 



(17) 



(Sari, Narayan & Piran 1996; Panaitescu & Kumar 2000; Sari & Esin 2001)^ is the Compton y- 
parameter of the plerion, which is the fractional energy gain of a photon when traveling through the 
plerion, due to upscattering by the relativistic electrons, and a = min[l, {jbrn / JbcY"'^] is the fraction 
of the internal energy in the electrons that is radiated away (Sari & Esin 2001). For our choice 
of parameters £{,5 <C ebe and there is fast cooling so that a = I. This implies Y^ ~ {fbe/^bs)^^'^ , 
and we shall use this relation in the following. The maximum LF is set by equating tc with the 
acceleration time, ~ 2'K^bmec/ qB^ (where q is the electric charge of the electron): 



The LF of an electron that cools on the adiabatic expansion time tgd (i-e. an electron for which 
tc = tsd) is given by 

= {l + t)Blart.. " «"''^273^i%^."ifi^rol53/?^i*s^d,o • (19) 

For < ti, where 

*1 ~ 1-1 V? ^'e,l/3^6S-3^rot.53/\-l , (^U) 

the cooling is so fast that Eq. (19) implies 7ftc < 1- This means that all electrons cool to non- 
relativistic random velocities within a time ~ 7bcisd and a distance ~ ^bcRb behind the termination 
shock of the pulsar wind. Of course, in this regime 7(,c < 1 no longer corresponds to the physical 
LF of the electrons^. We shall call this regime very fast cooling. 



''Sari and Esin pointed out that generally the factor a should be multiplied by the fluid velocity just behind the 
shock (in the shock frame), which for a relativistic shock with es <C 1, that is relevant for the pulsar wind termination 
shock, is /3 = 1/3. This factor of 1/3 should be divided by the ratio of the radiation flux and the photon energy 
density times c, which is 1/4 for an isotropic emission in the local rest frame of the emitting fluid. Together this gives 
a factor 4/3 which is close to 1, and is therefore neglected. 

^we formally obtain 76c < 1 because we used the approximation « 1 in the expression for the total synchrotron 
power of a single electron, under the assumption (that does not hold here) that the electrons are always relativistic. 
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In the case of a non-spherical (elongated) PWB, we can still use the same expressions, if we 
make the simple substitution described below. We define an effective radius, i?cfr = (SV/Att)^/'^ , so 
that a sphere of that radius will have the same volume, V, as the non-spherical PWB. Since the 
volume of the PWB determines its average number density and energy density (which determine the 
PWB emission), and since we expect the interior of the PWB to be roughly homogeneous (i.e. the 
local values are not very different from the mean values), then we can reproduce the expressions 
appropriate for a non-sphcrical model by simply replacing Rb with ReS- Since /3b appears in 
the various expressions only through Rb, we may achieve this task most easily by substituting 
f^b = 3i?cff/2ctsd everywhere. In order to illustrate the effects of a non-sphcrical geometry, we 
consider a particular example (which we also consider to be likely) of an elongated PWB with a 
polar radius much larger than the equatorial radius, Rp S> Req, where Rp/2 < R^s < Rp- If the 
surface mass density at the poles is sufficiently small, then the velocity of the SNR shell there can 
become close to c, so that Rp ^ ctgd- This would imply /?;, = 3Res/2ctsd ~ 1 



4.1. The Synchrotron Spectrum 



The characteristic synchrotron frequency of an electron with LF 75 is Ub = 7^f6o where Ubo 
SgBb/lQuieC. The synchrotron frequencies corresponding to 'jbrmlbM and 76c, respectively, are 



1/2 T-,-3/2 ^9/2 ,5/2 



l^bM 



5.1 X 10- (1 + ^)-S2%'4,l/3^6B:-3eJ/3^;ot:537^,4.5/36:TCd:0^ Hz 



.1/2 



1/2 



3-3/2-3/2 



(21) 



For tsd < ^2 [which is given in Eq. (35)] the synchrotron self absorption frequency Vbsa (which 
is treated below) is above the cooling frequency Vbci ^.nd the synchrotron flux density, F^,, peaks at 
Ubsa and consists of three power-law segments: 



^ i',max 



Vbc 



-1/2 



X < 



{v/vbsaf 
-1/2 



V < l^bsa 
I'bsa < fbm 



, {v\milvbsa) ^''^{vlvbm) < V < VbM 



For time separations between the SN and GRB events, t\ < tgd < ^3, where 

+ ^ 3/4 3/4 ^-l/2j^l/2 1/2 ^-3/2 

t3 ~ 65?72/3e(,;i/3e(,^^_34^_^^3£;,4^537j_4.5/3^__'i yr , 



(22) 



(23) 



we have 1 < ^bc < Ibm and the bubble is in the (moderately) fast cooling regime. For t2 < tgd < ^3) 
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Fi, peaks at Vc and is given by 



Fu{t2 < tsd < ts) 



= < 



saj 



V < Vbsa 
l'bc<V< ybm 



(24) 



For tsd > the bubble is in the slow cooling regime, where the spectrum peaks at Vbm and again 
consists of four power law segments: 

<Hsa 



{l^bsa/vbmY'^i'^/^b sal 
i^/l^bmf'^ 
{l^/l^bmf'-'^'^ 



^bsa < 1^ < l^bm 



(25) 



The peak synchrotron flux, F,^,max, for either fast or slow cooling, is 



F 



Ce-^rot-fi/jHiax (1 ~l~ -2-) 



42 (1 + ^)r?5;3'6^g_3ee,l/3<;':53Cl'^s"d:r7;;;4.5^Z28 ' (26) 



,3/2 ^-3/2,-3/2 -1 



where Pu,ma.x ~ i^csyn/^^syn, ^'e,syn = {4:/3)aTc{B'^ /8TT)-f'^ , V^yn = VbOll and di = 102^dL28 Cm is the 

the luminosity distance of the GRB. Synchrotron self-absorption (SSA) will cause a break in the 
spectrum at a frequency v^sa below which^ F^ oa v^. The absorption coefficient is given by: 



d 



Nile 



where 



is the spectral emissivity of a single electron. 



^syn/ 



for {v < Ui 



syn) 



(27) 



(28) 



The electron distribution, N{ie), is different for the fast cooling and slow cooling regimes. In 
the fast cooling case we have 

max(7ftc, 1) < 7e < 7bm 

(29) 

-2/ 



iV(7e) « A^o X < 



yibmi^ehi 



bm) 



-s-1 



Ibm < 7e < IbM 



®for PWBs, this applies also for the fast cooling regime, as opposed to the spectral slope of F^, oc v^^^^ (Granot, 
Piran & Sari 2000; Granot & Sari 2002) predicted for GRBs (both for the prompt emission and afterglow) since for 
the PWB the observer faces the back of the shell, i.e. the side farther from the shock, where the coolest electrons 
reside. 
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with A^o = Tftc^e- The local electron distribution depends on the distance behind the shock, I, 
and the distribution given above is averaged over I. Furthermore, it inchidcs only the relativistic 
electrons (i.e. for 7c < 1, most electrons would cool to non-relativistic random Lorentz factors, and 
are not included in this distribution). 

For u < max(t'6c, J^bo), all the electrons with 7e > max(7c, 1) contribute to the SSA and the 
absorption coefficient is 

17^ m™(7^^ ' , 1) . (30) 



a. 



If max(^';,c, i^w) < I' < i^bm we have 



and, finally, if > uiym then 



In the slow cooling case we have 



7e Ihm < 7e < Ihc 



iV(7e) « iVo X <( . (33) 

lbc{lehbc)~'~^ Ibc < 7e < 76M 

with No = {s — l)7b^^«e- The absorption coefficient for u < is given by 

"^-(3. + 2)8,^,.5/3,V3^^"^ (^4) 

The optical depth to SSA is given by Ti,^sa = Q;,yi?6 min (76c, 1), and i^^a is obtained by equating 
ru,sa = 1- For tsd < t2, where 

, ^-,0 3/75/21 4/21 ^2/21^11/21 -2/21^-31/21 

t2 ~ 12 r/2/3V,l/3^6B,-3^e,l/3^rot,53 7^,4.5 /^fo,-/ > (35) 

we have Vhsa > i^bc- Since typically ti < t2 < for igd < *2 we have either very fast cooling or 
moderately fast cooling and may use Eq. (31) to derive 

^ 1.9 X 10^° (1 + z)-\-'/'e-l{%e\,'l_,ej;i^^^^^^^ Hz . (36) 

For t2 < tsd < *3 we use Eq. (30) and obtain 

~ 1 1 V inl3 n -L ^^-1^1/2 6/5 1/2 7/10 ^3/5 ^9/5 -3/5 ^-24/5 -19/5 „ .07^, 
Z^6^a~l.lXlO [i + Z) a ' ?72/3Se,l/3^6B -3?e,l/3^rot,537tt,,4.5/'b-l *sd,0 ' (^'j 

Since the SSA frequency approaches Vbm only for very large time delays tsd 3> ^3, we do not consider 
the case Vbsa > ^bm- In the slow cooling regime (tgd > ^3) we have 

7/ ~ 9 Q V in9 n -U r^-lr.-'^/S^^/S ;t8/5 p4/5 ^"8/5^-9/5 -9/5 .oo\ 
Ui,sa ~ ^-9 X 10 (1 + Z) ?72/3 ^6B,-3^e,l/3^rot,537^,4.5P6,-l ^sd,0 ' (^O) 
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4.2. The SSC Spectrum 



The SSC spectrum has a similar shape to the synchrotron spectrum, and we approximate it 
as comprising of broken power-laws with characteristic frequencies. For the different ranges in tsd, 
it assumes the following forms: 



,SC* 
bsa 



(,/,5C)(2-.)/2 



<^< ^bm 



(39) 



I' > Ui 



uF^C{t2 < tsd < ^3) 



Yb ^bmPi 



= < 



where = max(l, 7^^)fbsa, and 



(,/,^^C)(3-.)/2 
(,/,5C)(2-.)/2 



sc 

bm 



SC 
bsa 



^bs^a < ^ < ^bc 



^bc < ^ < ^fm 



(40) 



sc 

bm 



sc 

bsa 



^b^a <^< ^bm 



^bm <^< ^bc 



(41) 



sc 

be 



SC 



^In, ^ 2.1 X 10^ (1 + e2^^^/3e,^^;:3i?,,;;5l,/3,7i^t;,7 Hz 



3/2 p-7/2 ^21/2 ,13/2 



SC 



iL^bm « 6.4 X 10^" (1 + ^)-'^2/34,l/344'-3^e;i%^rot:537^,4.5Crisd:o' Hz , (42) 



1/2 4 



-3/2-3/2 



and 



/'^2„ ^ 1.4xl0i° „-7/6^-2 ,-7/6 ,-5/6 ^1/3 ^-5/3 -l/3ol6/3,3 „ . . y. , 

Tbc'^b^a ~ ~(T+ir ^2/3Se,l/3^6B -3?e,l/3^rot,537^«,4.5P6-l*sd,0 *1 < < *2 



^bsa 



„2,, 8.0x1012 ^-1/2^-4/5,-1/5 ,-3/10^3/5 p-1/5 -3/5^6/5 1/5 „ 

" - " ''2/3 ^fee, 1/3^6-8, -3^6, l/3^rot,53'ui,4.5^\-l''sd,0 ''2 ^ (-sd ^ 



Tbc'^fesO ~ (1 + 2) 



L Ibm^bsa ~ (1+2) '72/3^66,1/3^65 _3?e,l/3^rot,537^,4.5P6,-l *sd,0 *sd > 



1/5 f-2/5p4/5 2/5 ^-9/5,-9/5 



(43) 



The peak of vF^^ is simply times the peak of the synchrotron vFy, which for fast cooling is 
given by 



J^fem^i/fc^ - Z^6m^i/,maxy^~ 5.1 X 10 ^^a ^^^%/3Se,V3^&B%-^''°t'53*sdVL28 ' l^^) 

and for slow cooling is given by this expression multiplied by the factor a. 
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4.3. The Klein-Nishina Effect 



When the energy of the seed (synchrotron) photon in the rest frame of the scattering electron 
exceeds the electrons rest energy, jehi^'"^^"" > meC^, then we move from the Thomson limit to the 
Klein-Nishina (KN) regime, where there is a reduction in the scattering cross section and the 
electron recoil becomes important. The corresponding frequencies where the KN limit is reached 
for the seed synchrotron photon and upscattered (SC) photon are given by 



sc 

'^KN_ _ syn 
— lei^KN 

le 



i^KNfi = ^ = ^ = 1-23 X lO^o Hz , 
h Ac 



(45) 



where Ac is the Compton wavelength. As the energy increase of the photon due to the scattering 
cannot exceed that of the electron, the photon energy of hvj^^ = 7emeC^ sets the natural upper 
limit to the frequency of photons that are upscattered by electrons with a LF 7e. Typically, either 
the fractional energy gain in successive scatterings is Y ^Tll <S 1, or the KN limit is reached for 
the second scattering, 'y^hv^'^'^ > mgC^, thus allowing us to ignore multiple scatterings. 

In KN effects typically become important only at i^'^'^ > max{i'^^,u^^) in the fast cooling 
regime, or at > in the slow cooling regime. We therefor restrict our discussion to these 
frequency ranges of the SSC spectrum. Within this frequency range, the SSC flux density at a 
given frequency consists of roughly equal contributions from seed synchrotron photons that extend 
over a finite range of frequencies, that are upscattered by electrons within a finite range of LFs 75, 
so that 7^2^*^"" = u^'^' = const. For this reason, a significant change in the SSC flux at a certain 
frequency u'^^ , due to KN effects, will occur only when all the electrons that contribute significantly 
to this frequency reach the KN limit. Since i^xifi'^b) 7;,, this occurs when the electron with the 
maximal 75 that still contributes significantly to u^'-^, Jmaxii^^'''')-, reaches the KN limit. For fast 
cooling 



'Jmax (^) — * 



^ > IbM^bm 



(46) 



^, IbM 

where viyp = max(i/65a, i/{,c), while for slow cooling 



76c 



^vlvbc < ^ < ^bM^bc 



(47) 



IbM 



^ > IbM^bc 
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where uj^ga < Vhm- The KN effects become important at tlie frequency ffj^i that satisfies^ v 
lmax{i^)^KN,Q- For fast cooling 



^kNA — \ 



(48) 



while for slow cooling 



2 ; 

< 7bMl^KN,0 

Ibc^KNfl 
K IbMl^KNfi 



^b^ < 4na < IbMnm 



4n,1 > IbM^bm 



^b^ < 4n,1 < Ibcnm 



sc 



Xc^bm < ^kNA < 



,sc 

^bc 



(49) 



4c < 4n,1 < IbM^bc 



4n,1 > '^bM^bc 



At < v^Ni effects are unimportant and the SSC spectrum is given by equations (39-41). 
At 1/ > f^Ni effects become important. In order to calculate in this range we first need 
to estimate the scattering optical depth re(76) of electrons with LF 7e ~ 76, which is of the same 
order as the total optical depth of electrons with 7e > 7^: 



rlbM 

•lib 



(50) 



where Nipi^) is given in equations (29) and (33) for fast cooling and for slow cooling, respectively. 
For fast cooling we obtain 



Ibche 



Teile) ~ Tt X < 



. Ibcliyfn 7e 



while for slow cooling 



' {lehbm 



A~s 



Teile) TT X 1 



. Ibcliyfji 7e 



76c < 7e < 76m 
76m < 7e < 76M 

76m < 7e < 76c 
76c < 7e < 76M 



(51) 



(52) 



''from equations (46) and (47) it is evident that there is always exactly one solution to this equation, while the 
explicit form of this solution depends on the frequency range in which it is obtained. 
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Kiv > ly^^i the SSC flux density is dominated by the contribution from electrons with i^xNilb) ~ ^ 
and is given by 

F,'''^Te{je = l^/l^KN,0)F^r ■ (53) 

The frequency dependence arising from the first term on the r.h.s of equation (53) is just the 7e 
dependence of re(7e) [that appears explicitly in equations (51) and (52)], while the second term 
introduces a frequency dependence of where F,^^" oc at v^^^ = '^xato/^^'^- 

In all cases the largest frequency that the SSC spectrum reaches is min[z/^.^(7{,M)5 7feM^fe-^]- 
For T-'xiii^hM) > IhM^bM the SSC spectrum is given by equations (39-41), with no changes. For 
7^jy^ max ( i/ftc) < T^Kki'ybM) < JbM^bM, the SSC spectrum is the same as in equations (39-41) 
up to fxSiijbM) where it sharply ends [note that v^fj^i'^bM) = ^ in this region]. 

For fast cooling with Jbm'^bp < ^^Iw,! < IbM^bm we have F^^ oc v^f^^^ for v^fj^ i < v < 
^xki^bM), and the spectrum ends at ^'^*^(76m)- Immediately above z^^ati we have (3 = —1/2 
while /? might change to 1/3 (for Ubc > Vbsa) or 2 (for Ubc < Vbsa) at a higher frequency, i^xNo/i^bpi 
producing a spectral break at this frequency, if i^xj^^o/i'bp < ^RNhbu) [i-e. if llu^bp > '^KhilbM)]- 
For fast cooling with 7^^1/5^ < zyf'^-^ < -fl^.^bp wc have F^^ oc u''^-'^ for i^ff^^i <v < u^%j{-1bm), 
and F^^ oc u~^^'^ for i^ff^i-Jb-m) < v < ^liiiilbM), where for Vbsa > I'bc we have (3 = 2, while for 
^bc > i^bsa we have (3 = 1/3 immediately above i^^k 1^ which may change to /3 = 2 at o/i^bsa if 

For slow cooling with ^l^Mbm < ^icN,i < IbM^bc wc have F^^ oc v^^^^ for v^k,! < ^ < 
^xifi'ybM), and the spectrum ends at i^xki'^bM)- Immediately above i^xk i have /3 = (1 — s)/2 
while j3 might change to 1/3 at a higher frequency, i^l^jsj Q/i^bm, producing a spectral break at this 
frequency, if VKN.o/^bm < t^KkilbM)- For slow cooling with 7^/^ < vffj^ i < ^Ic^bm we have F^'^ oc 
j^i-s-p foj, < J, < i^Kkhbc), and F^^ oc i^-'-^ for i^fg^i-fbc) < < t^Kkilbu), where /? = 1/3 

immediately above i^^N i, and may change to /3 = 2 at f^Nfi/^bsa if ^KN,o/^bsa < J^|w(76m)- 



4.4. Opacity to Pair Production 

High energy photons emitted in the PWB, either by the plerion itself or by the prompt GRB 
or afterglow that occur inside the PWB, may interact with lower energy photons of the strong 
radiation field of the plerion to create pairs. For sufficiently high photon energies, the optical 
depth to this process, r-y-y, may exceed unity, so that they could not escape and reach the observer. 
We now calculate the photon energy e (in units of mgC^) for which r-y^(e) = 1. This sets the 
maximal photon energy that will not be effected by this process. 

The radiation field of the plerion is roughly homogeneous and isotropic within the largest 

radius where the radiation is emitted, which we parameterize as fRb- For a fast cooling PWB 
{tsd ^ ^3 ~ 65 yr), the radiation is emitted within a thin layer behind the wind termination shock, 
at the radius Rs, so that / = Rg/Rb For an adiabatic bubble, as the one considered here, the 
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value of this ratio ranges between 0.2 and 0.5 (KG). For a slow cooling PWB {tgd > ^ yr) 
the radiation is emitted from the whole volume of the PWB, and / = 1. The internal shocks 
that produce the prompt GRB emission take place at a radius smaller than fRb, and the energy 
density of the plerion radiation field is Uph ~ a^e-£'rot/*sd27r(/i?ft)^c. However, the relevant target 
photons for pair production with high energy photons are the synchrotron photons, since the 
synchrotron component is dominant at low energies. Therefore, we should use the energy density 
of the synchrotron photons, Usyn = Uph/{1 + 1ft). As we shall see in §6.3, the afterglow emission 
typically occurs at i? > fRb, where the radiation field is not homogeneous, but rather drops as 
Uph oc R^^, and is not isotropic, causing a smaller typical angle between the trajectories of the two 
photons that could possibly produce an pair. Both effects reduce r-y-y for the afterglow emission, 
compared to that for the prompt GRB or the plerion emission itself, that we calculate below. 

The number density of synchrotron photons, n^, per unit dimensionless photon energy e, may 
be obtained from the shape of the synchrotron spectrum, and its normalization is 

O^e-E'rot 2 /" i 2 2 /r;j\ 

ct^2^{fR,ni + Y,) = ^^y- = "^^^ J = ^"^^^ ' 

where Uph is the photon energy density of the plerion and g = 2(s — l)/(s — 2). The optical depth 
to pair production is given by 

I \ -PD -1 ^A'y'l 2)/2 s/2 (r-r-\ 

r-y^(e) « axfRbe Wi/e = aTjRb —^m > (55) 



and r^-y(e) = 1 is satisfied for 

_ eruec _ 2.4 GeV / A/3\ " -5/22 3/11 -1/2 ^-8/11^-21/22 -2/11 .23/22 43/22 . . 

nv^^ - ^-|^ _|_ ~ ^-^ _|_ V / ''^ ^&e,l/3^6S,-3?e,l/3 -^rot,53 Iw,4.5 Pb-1 hdfi ' K^^' 

where /1/3 = //(1/3). It can be seen from equation (56), that opacity to pair production becomes 
important only at very high photon energies, and is larger for smaller tgd- 

For an elongated PWB, the radiation field can generally have a rather different structure, 
resulting in a different expression for v^j. However, one can imagine a simple scenario where the 
structure of the radiation field is similar to the spherical case, and equation (56) is still applicable 
with simple substitutions. This can occur if the PWB is fast cooling and most of the radiation 

is emitted just behind the wind termination shock, and the latter is roughly spherical, with a 
radius, Rg, similar to the equatorial radius, i?cq- In this case we should make the usual substitution 
Pb = 3i?efr/2ctsd5 as well as / = Rg/ReS Req/ReS- For example, with w 1 and Req/ReS ~ 0.1 
this would increase u^j by a factor of ~ 3. 



4.5. Prospects for Direct Detection 



An important prediction of this model is a strong radiation field within the PWB. We now 
examine the possibility of directly observing the radiation emanating from the PWB during the 
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time between the SN and the GRB events. For time separations tgd smaher than 



0.4 (^] 



1/2 



(57) 



the SNR shell has a Thomson optical depth larger than unity, and would therefore obscure the 
radiation emitted within the PWB. For tgd < U the emission due to the radioactive decay of Ni 
and Co in the SNR shell, might be observed, as in a regular supernova. However, even at the peak 
of the supernova emission, it will be hard to detect at a cosmological distance. This difficulty is also 

present in ongoing searches for high redshift supernovae, where like in our case, random patches of 
the sky need to be searched, as there is no prompt GRB and afterglow emission to tell us where 
and when to look. 

If there is considerable clumping of the SNR shell before this time, then equation (57) gives 
the time when the average optical depth of the SNR equals 1, while for regions of the shell with a 
less than average density the optical depth can drop below unity at a somewhat earlier time. This 
constraint may also be eased if the geometry of the PWB is not spherical (an elongated PWB), 
and the mean density of the SNR shell is significantly smaller near the poles compared to near the 
equator, and our line of sight is near one of the poles (as is required in order to see the prompt 
7-ray emission from a jetted GRB). 

Once the SNR shell becomes optically thin, the PWB radiation may be detected if the flux 
that arrives at the observer is sufficiently large. For simplicity, we calculate the flux at the time 
of the GRB explosion, tg^ after the supernova event, since the relevant quantities scale as power 
laws with the time t after the supernova, and the system spends most of its (logarithmic) time 
near t = tgd- For concreteness, we consider the observed flux at the radio, optical and X-ray, for 
a typical frequency in each of these frequency ranges: frad = 5 x 10^ Hz, i/op = 5 x 10^^ Hz and 
z/x = 10^^ Hz. 



For the radio we have frad < Hsa for tgd < ^sa.rad, I'bsa < Z^rad < Z^6c for *so,rad < tsd < ic,rad 

and < I'rad < ^bm for isd > *c,rad, where 



(-sa.rad — o-^ i" Z) tfee,l/3*'6-B ,-3'?e,l/3-^rot,53 /w),4.5/^6 -1 '^9.7 



t A - 20(1 + 2)2/5„3/5 2/5 1/5 ^3/5 ^"9/5 2/5 
fc,rad — ^U[L-\-Zj ^2/3^6,1/3^65 -3-^rot,53/^6,-l ^14.7 ' 



p3/5 



-9/5 2/5 



(58) 
(59) 



and 



psyn ^ ) 



34 n3y -1/4 -1/12 -1/6 ^1/6 ^-1/12 -1/6^29/12 9/4 o ^ 
(l+2)-3 '/2/3 ^6e,l/3^6B ,-3?e,l/3-^rot,53 7«,,4.5/^6-l ^sd,o'^9.7" 



2 

L28 



(l+2)-V2 '?2/3 Se,l/3^6B,-3^e,l/3-'^rot,53'Ttt,,4.5/^6,-l''sd,0 ^^9.7 ^L28 '^sa,rad < ^sd < l^c,rad 



-1/4 -1/2 1/4 



-,3/4 



1 «3/4 ,-1/4, -1/2 ,_2 



tsd < t. 



sa.rad 



0.50 Jy 1/3 2/3 ^ ^2 -1 /3-3 ..-7/3 1/3,-2 . . . 

(l+5p473 ^2/3^^6,1/3^65 _3?e,l/3-C'rot,537M,,4.5P6,-l''sd,0 ^9.7 "L28 ^sd > Ic,md 



3 ,-7/3, 1/3 ,-2 



(60) 
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For tsd < tm,op we have t'op < ^bmi while for tg(j > tj^^op, the ordering is reversed Ugp > t^brm 
where t^,op is given by 

f - 22 fl + 2r2/3r7^/3 4/3 1/3 ^-4/3p.l/3 4/3 ^-l -2/3 .... 
Cm,op — U.ZZ (^i -I- ^2/3^6e,l/3^6B-3^e,l/3-^rot,537M,,4.5n -1^14.7 ■ 

The optical flux is dominated by synchrotron emission and is given by 

, _3JjxJy -1/4 -1/2 1/4 ^ ^3/4 l ^3/4 -1/4 -1/2,-2 , 

(H-z)-V2 '/2/3 ^6e,l/3^6S -3?e,l/3-^rot,537u,,4.5P6,-l'^sd,0 '^14.7 "L28 ''sd ^ ^m.op 

psyn _ J 

I 0.82 MJy 5/4 7/10 11/20 ^-1/5 p21/20 1/5 ^-3/20 -23/20 -11/10 ,-2 + ^ + 

^ (l+z)i/io '/2/3^be,l/3*'6-B,-3?e,l/3-^rot,537t(;,4.5'^b,-l ''sd,0 '^14.7 "L28 ''sd ^ ''"i.oP 

(62) 

The X-ray is typically above Uf^m, and for the synchrotron emission we have 

,,Z7^?/n _ X 10~^^ ^5/4,7/10 11/20 ^-1/5 ^21/20 1/5 ^-3/20 -23/20 -1/10 ,-2 erg 

'^-^i^ ~ _|_ ^^1/10 ''2/3^6e,l/3^6B,-3Se,l/3-^rot,537u,,4.5P6,-l '^sd,0 ^^18 ^28 ^^^2 g " \^'^) 

For sufficiently small values igd we have < i^x < As tgd increases, then decreases 

below the X-ray for > tf^y^-, and i/^'^ grows above the X-ray for t^A > ifx' where the relative 
ordering of tf^y^ and if ^ depends on the values of the other parameters. For our fiducial values we 
have 

*m,X - ^0(^1 -l-Zj '72/3^6e,l/3^6m^6B,-3Se,l/3-^rot,537tt,,4.5M6,-1^18 ' I04j 

fSC _ en f1 + ?12/13 4/13 7/13 4/13 3/13 ^7/13 ^-21/13 2/13 ..^.^ 

For this ordering of these two times we have z/x > ™ax(f^, i/^'^) for t^*-^ < igj < tf^, and 
< < t'fc^ for tsd > if,x- For the other ordering, t^'^ < t^^, we have ux < min(i/f^, u^^) 
for < tsd < t^x '^bm < ^ < '^bF for tsd > t^x- The X-ray is always below the KN limit, 
and the SSC vFi^, is given by 



vF 

ly J. 



sc 



(66) 



. 3.7-10-^-^ ^-5/4^-1 -1/4^2 p3/4 2 ^7/4 -1/4 1/2,-2 , /.5C^5C\ 

(l+2)-i/2^2/3 ^6e,l/3^6B -3Se,l/3-^rot,537u,,4.5P6,-l'^sd,0 ^^18 "L28 '^sd ^ min(,rc,X' ^^m.xJ 

1.8-10-13 ^29/20^7/5 1/20 ,-2/5 ^21/20 2/5 ^17/20 -23/20 -1/10 ,-2 +SC ^ + ^4. 

(1+2)1/10 '/2/3 *^be,l/3^6B,-3?e,l/3^rot,53 /«),4.5/^6,-l ''sd.O ^18 "L28 '^m,X ''sd ^3 

2.8-10-1-1 „ 16/15 -11/180^-8/45 p5/6 8/45 ^31/60 -127/180 -1/10 ,-2 . . .SC 

(1+2)1/10 '/2/3e6e,l/3^&B -3 ^£,1/3 -^rot,537«,,4.5A^6,-l '^sd,0 ^^18 "L28 '^3 ^ ''sd *^ ^c.X 

2.0-10-5 ^5/3,2/3 ^2 p.11/3 2 ^-7.-17/3 4/3,-2 ^SC ^ f ^ f 

(l+2)-4/3 '/2/3^6e,l/3^*''S,-3?e,l/3'^rot,537«,,4.5P6-l'^sd,0 ^^18 "L28 ''c.X ^ '^sd ^ C3 

3.010-1° „-l -4/3 1/3 ,10/3 p7/3 ^-10/3^-4 .-3 4/3,-2 U < f , <r fSC 

(l+z)-*/3 ''2/3^6e,l/3 6B,-3?e,l/3 rot,53 7«,,4.5 /-^fe,-l''sd,o'^18 "L28 ''S ^ ^-sd ^ ^m^x 

1.2-10-^ „16/5 12/5 4/5 ^-2/5 ^14/5 2/5 ^-22/5 -22/5 2/5 ,-2 f . ~> mf^^cit^^^ t^^^ \ 

(l+z)-2/6 ''2/3 ^6e,l/3^6B,-3Se,l/3-^rot,537«,,4.5P6,-l ^sd,0 '^18"L28 ^sd ^ "laxi^C^ x, 

The emission from the PWB calculated above is at the time of the GRB explosion, tgd- How- 
ever, an important question one needs to address is for how long after the onset of the GRB 
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explosion will the radiation from the plerion persist. Once the SMNS collapses to a black hole, 
giving rise to the GRB event, the pulsar wind stops abruptly, and no new electrons, freshly acceler- 
ated at the wind termination shock, are injected into the plerion from this point on. The electrons 
in the PWB begin to cool radiatively, and adiabatic cooling of the electrons become important on 
the dynamical time of the plerion, tg^ (where this is a good estimate of the dynamical time also 
for an elongated PWB). Once the last accelerated electrons cool below the Lorentz factor at which 
they radiate at some observed frequency no more radiation is emitted at that frequency. This 
cooling time is given by 

+ - a Q n -U ^^-1/2^-3/4 -1/2 -1/4 ^9/4 ^-3/4 9/4 -1/2 , , . 

tcool - 8-9 (1 + ' '?2/3 ^fee,l/3^bi?-3/'b,-l-^rot,53^sd,0'^9.7 ^lays , (b7j 

where the above numerical coefficient is for the radio, while for the optical and X-ray the numerical 
coefficient is 40 min and 24 s, respectively. For our simple example of an elongated PWB, where 
/3b « 1, the cooling time is about two orders of magnitude larger. Due to the strong dependence 
on tsd, icooi can become quite large for large values of tgd; especially for radio frequencies. 

A possibly more severe constraint arises from the geometrical time delay in the arrival of 
photons to the observer, from the different parts of the PWB 

ig ~ = 24.2 /36,-itsd,o days . (68) 

This is to say, that even if the emission from the PWB would stop at once, with the onset of the 
GRB, the radiation would still reach the observer for after the GRB, just due to the geometrical 
time delay in the arrival of photons to the observer from the far side of the PWB, compared to the 
side facing the observer. For t^d < ^3, where t-j, is typically rather large [see equation (23)], the PWB 
is in the fast cooling regime, and most of the emission occurs near the radius of the termination 
shock, Rg. This reduces tg by a factor of Rg/Rb- However, this is not expected to account for 
more than a factor of ~ 5 (KG). For an elongated PWB we have tg ~ Rg/c for tgd ^ ^3) while for 
isd > ^3 we have tg Rp/c ^ tgd- The emission from the PWB should persist for an observed time 
of tpi ~ (1 -|- z) max[tg, min(tsdi ^cool)] after the GRB, which as can be seen from equations (67) 
and (68) should be at least a few days after the GRB, but can also be much larger {tpi > tg, see 
equation 68). 

The plerion emission at the radio, optical and X-ray bands is shown in Figure 1. For reference, 
we also show the times (below which the Thomson optical depth is larger than unity), tp^ (below 
which iron line features can appear in the X-ray spectrum of the afterglow) and tisM (for which 
the effective density of the PWB is similar to that of a typical ISM, i.e. 1 cm"^). In the radio, the 
typical limiting flux for detection is ~ 0.1 mJy, and upper limits at this flux level, at a time t < tpi 
after the GRB, would exclude 2 < tg^.o ^ 20 for ^ 10^, while for 2 < tsd,o ^ 65 this would imply 
Iw ^ 10^. Values of tgd ^ 2 yr or tgd ^ 65 yr are hard to constrain with the radio. 

Optical upper limits at the R-band at the level of 24 — 25th magnitude {F,y < 0.5/xJy) would 
imply jyj > 10^ for tgd < tpe, i-e. for afterglows with iron lines. More stringent upper limits, may 
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provide more severe constraints. For example, an upper limit of 0.01 /xJy (i.e. R=28.6 which might 
be reached with HST) implies tsd ^ iiSM- 

In the X-ray, for 7^ < 10^-^ the SSC emission dominates over the synchrotron emission for all 
tsd ^ tr, and for j^j = 10^ it dominates for tsd ^ 7 yr. The typical limiting flux for detection in 
the X-ray is a few xlO~^^ erg cm~^ s~^, and upper limits at this level may imply > 10^ for 
tsd ^ ^Fe- For tsd ^ tpe, such upper limits cannot provide any useful constraints. 



5. Effects on the Prompt GRB Emission 

The prompt gamma-ray emission is believed to arise from internal shocks within the GRB 
outflow, due to variability in its Lorentz factor T (Rees Sz Meszaros 1994; Sari Sz Piran 1997). In 
order for this process to be efficient, it needs to occur before the ejecta is significantly decelerated 
by the ambient medium. Therefore, the main effect that a plerionic environment may have on the 
prompt GRB stage is through inverse Compton scattering of the photons from the external plerion 
radiation field (which we shall refer to as external Compton, or EC). 

The external (plerion) radiation field in the local rest frame of the shocked shells is [7^^ = 
V^Uph^ext where Uph,ext is given by equation (54). The electrons radiativcly cool by the combination 
of the synchrotron, SSC and EC processes the timcscales for which are (in the comoving frame) 
igyn ~ 67rmec/(7T-B'^7, t'g^ = t'syn/Y and t'^^ = t'sy^/X, where Y is the Compton y-parameter and 

U' 

is the ratio of the energy density of the external radiation field and the magnetic field in the 
local rest frame of the shocked shells, which is also the ratio between the energies in the EC and 
synchrotron components. As can be seen from equation (69), X <C 1 for typical parameters. In 
order to have X > 1, we need tsd i$ 1 yr and F > 10^. For an elongated PWB of the type described 
just before §4.1, with / ~ Rg/ReS ~ 0.1 and Pf, ^ 1, X is roughly the same. 

The EC component is due to the scattering of external photons from the plerion radiation 
field by the electrons from the GRB ejecta. This scattering can be done either by hot (relativistic) 
electrons, or by cold (non-relativistic) electrons, the latter being either in cold shells or cold portions 
of colliding shells (in either regions before the shock, or at a distance larger than 7cA' behind the 
shock for 7c < 1). 

For the colliding shells, we assume that Thomson optical depth, tt, is smaller than 1. We 
provide detailed expressions for one representative plerion spectrum, the one given in equation 
(22), that is relevant for 4d < ^2 ~ 12 yr (see equation 35), which is of most interest. Similar 
expressions for the other plerion spectra can be readily derived in a similar manner. The plerion 
SSC emission in this regime has a peak for the uFi^ at > 10^^ Hz for tgd ^ 1 yr, and will 
therefore be above the KN cutoff for both hot and cold electrons in the outflowing shells, and its 
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contribution to the EC emission can be neglected. The resulting EC spectrum due to scattering 
by the hot electrons is 



pEC 

■ Xxl 



{-l-l'^fl^ -fa^<-<-^^ ■ (70) 



where = max (7^, 1)1^6^0 and 

EC _ t2 2 ..^ 1 n in23 5/2 2 1/2 >-2 i-,l/2 2 ^-3/2 -3/2 2-p2 tt 

— ^ Im^bm ^ i-->i X i^'J V2/3^'be,l/3^'bB-3^e,l/3^Tot,53^w,4.5Pb,~-l'^sdfl^e^2.5^'^^ 

SC _ 2 / ^Qn^in27 -5/2 _2 -1/2 ^2 171-I/2 -2 o-3/2,3/2 
^^JCAT,! - ^KNfl/^bm ~ -3-^ X iU ??2/3 ^6e,l/3^6B,-3^e,l/3^rot,537«,vl.5P6 -1 f^sd.O ' 

^Wn{1m) = rjMl^KN,0 ^ 6.6 X 10^^(1 +Y + X)-V2ey4,^l/4^-V4r5/2^1/2^ _ ^^^^ 

If < ^kn{'^m) then we have z/F,, oc 1^^2-5 f^j. ^,|C|^^^ < < i'kn{'^m)- The peak of the vF^'^ 

spectrum, of the EC component from hot electrons is, is ~ 10~^°(X/10~^) ergcm^^s"-*^ and is a 
factor of X ~ 10~^ (for t^d ~ 1 yr) smaller than that of the synchrotron component, and therefore 
might be detected only for extreme parameters (isd,o i$ 1) T 10^). 

For the scattering by cold electrons, the optical depth is approximately 

TT = at<A' « Qme-^L^2V2lt-^_2 , (72) 

and the Compton y-parameter is 1" = min(rT, 1)F^. This implies that all frequencies of the plerion 
spectrum are shifted upwards by a factor of F^, and the corresponding flux density, Fi,, should be 
multiplied by min(rT, 1). The only exception to this simple re-scaling is that the spectral slope 
below the upscattered self absorption frequency will be Fi^ oc v, instead of the in the plerion 
spectrum. The peak of the uF^^ spectrum will be at 

T'u,^ ^ 5.1 X 10^« (1 + ^)-\2/34,l/34!-3^e":i%^r5537^,4.5/36:!('isto' , (73) 

which is typically at the hard X-ray or soft gamma-ray for tgd ^ 1 yr. However, the peak of uF^^ 
for this component is a factor of F < F^ 10^ larger than ui^^F^^^ < 10~^^ ergcm^^s"-*^ and is 
therefore < 10~^ ergcm~^ s~^, and will typically hide below the standard GRB emission. 

Another possible effect of the plerion radiation field is that photons with energy > 1 q GeV 
cannot escape the emission region due to a large opacity to pair production (r^^ > 1). Therefore, 
all the components of the prompt GRB emission, including synchrotron, SSC and EC, will have an 
upper cutoff at this photon energy. 

Finally, we consider the effect of the Compton drag due to the plerion radiation field on 
the GRB outflow^. The effect of Compton drag in GRBs was considered in the context of the 



*The emission from the initial supernova itself always contributes much less to the total radiation field inside the 
PWB, and may therefore be neglected. 
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collapsar model, where the radiation comes from the walls of a funnel along the rotational axis of 
the progenitor star (Ghisellini et al. 2000; Lazzati ct al. 2000). The rate of energy loss of a shell 
of initial Lorentz factor Fq rest mass M and solid angle Qj , is given by 



H ~ ~dt 



Mc^ = -njR^cr^Uphmm{TT, 1) 



(74) 



where t ~ R/c is the lab frame time, tt = {R/Rr)~'^ is the Thomson optical depth of the shell, 
Rr is the radius where this optical depth drops below unity^. We render Eq. 74 dimensionless by 
introducing R = R/ Rt = Tr^ 



-1/2 



This gives 



where 



dim 

dR 



T{R) 



^min(^^l) 



A = 



(1 + AToR^/3)-^ R < 1 



l + Aro{R-2/S) 



-1 



(75) 



(76) 



R>1 



Rr = 3.3 X 10^2 ^;l/2^1/2r^ 1/2^1/2^ cm, 



(77) 



where Sj is the fraction of the kinetic luminosity of the GRB outflow that is converted into the 
gamma-ray emission. As can be seen from Eqs. (76) and (77), ATq ^ 1 and therefore T{Rt) = Fq, 
while for R ^ R^ the fractional change in T is given by ATqR. If the radius, -Red; &t which 
deceleration due to Compton drag becomes significant, is larger than the deceleration radius, i?dcc ~ 
fRb, due to the sweeping up ofthe PWB material, then the deceleration due to Compton drag is at 
most comparable (and never dominant) to the deceleration due to the ambient medium, for^*^ A; < 1. 
Therefore, Compton drag will have a significant effect on the decelleration only if Red < Rdec, which 
for our fiducial values may be expressed as tgd < 0.3/^/^ yr. For such low values of tgd the SNR 
shell is still optically thick to Compton scattering (tsd < tr), so that we do not expect to see the 
GRB or afterglow emission. We conclude that the deceleration of the GRB ejecta due to Compton 
drag is negligible for relevant values of igd- 



We have used the total photon energy density ofthe PWB, Uph, that includes the SSC component, even though 
for tsd ^ 30 yr, most of the energy in the SSC component is in photons that are above the Klein-Nishina limit, and 
would therefore have a reduced cross section for scattering. Since we show that Compton drag is unimportant even 
without taking into account the reduced cross section, this effect can only strengthen our conclusion. 

^"This is since the Lorentz factor decreases with radius as due to Compton drag and as J?(*'~^)/^ due to the 
ambient medium. 
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6. Effects on the Afterglow Emission 



At a time igd after the supernova event, the SMNS cohapses and triggers the GRB explosion, 
sending a fireball and relativistic blast wave into the PWB. When the GRB ejecta has swept up 
enough of the outlying material, it is decelerated, and it drives a strong relativistic shock into the 
external medium, that is responsible for the afterglow (AG) emission. 

In this section we study the observational consequences of the plerionic environment inside 
the PWB, that are different from the standard "cold", weakly magnetized proton-electron external 
medium. One of the advantages of having the PWB as the environment for the GRB afterglow is 
that it naturally yields high values of eg and (the fraction of the internal energy in the electrons 
and in the magnetic field, respectively) behind the AG shock (KG). High values of eg are expected 
from the fact that relativistic pulsar-type winds are likely dominated by an electron-positron com- 
ponent, whereas significant values of should naturally occur if the winds are characterized by 
a high magnetization parameter. We expect ^ e^B, and use the same fiducial value for these 
two parameters. The electrons in the PWB are typically colder than the protons by the time the 
afterglow shock arrives, and most of the energy is in the internal energy of the hot protons. This 
might suggest that eg can be slightly smaller than e^g and motivates us to use eg = eg,_i/10 for our 
fiducial values. 

The values of the physical quantities behind the AG shock can be determined from the appro- 
priate generalizations of the hydrodynamic conditions used in the case of a "cold" medium taking 
into account the fact that the preshock gas is now "hot" and should be well described by a rel- 
ativistic equation of state, p = e/3 = w/A, where w is the enthalpy density and p is the particle 
pressure. In the following we largely follow the analysis of KG . The deceleration of the AG shock is 
determined by the total enthalpy of the external medium, wtoti which includes contributions from 
the particles and the magnetic field enthalpy, B'^/Att, where the latter contribution is negligible for 
our choice of parameters {cb ^ !)■ This make it convenient to define an "equivalent" hydrogen 
number density nH,equiv = Wiox/mpC? ^ w/nipCp = (4/3)e/mpC^, in analogy with the traditional 
parameterization of the external medium enthalpy density, w = numpCp, that is relevant for a 
standard ISM or stellar- wind environment. 

In general both the energy and the electron number density may be function of the distance r 
from the center of the PWB and can be parametrized as 



e(r) 



Aer 



(3 - K)r]E^ot 



(78) 



ne(r) = AnV 



-k 



(3 - k)Ne 



(79) 



where iVg = Ngtsd is the total number of the electrons in the PWB and iVg is given in Eq. 10. 
When a large fraction of the energy density in the PWB goes to the proton component we have 
rj ^ 1 and we can expect both nn.equiv and Ue to have a similar radial dependence, i.e. k = k*. 
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The expected values of k typically ranges between k = 0, similar the the ISM, and k = 1, that 
is intermediate between an ISM and a stellar wind (KG). For an elongated PWB, things can get 
much more complicated, since e and rig can depend not only on r but also on the angle 9 from the 
polar axis. However, if the 6 dependence within the opening angle of the GRB jet is small, and the 
dependence on r may be reasonably approximated by a power law, then our formalism still holds 
for A; 0, with the usual substitution of Pf, = 3Res/2ctsd (see the beginning of §4). For A; > we 
also need to change the normalization in equation (78) and (79) accordingly. 

The expressions for the radius Rag and the Lorentz factor Fag of the shocked AG material 
can be derived using energy conservation 



i-Rag 

AG / nH,equiv?^^pC^47rr^dr , 
JRs 



where E\so is the isotropic equivalent energy of the AG shock, and using the relation 

Rag 



t 



(80) 



(81) 



where t is the observed time. We obtain 



Rag = 



^EisoCtR'i^ 



3-fc 



r]E, 



rot 



l/(4-fe) 



Rag 
Act 



(fc-3) 



= 



k 



(82) 



(83) 



The postshock energy and particle density (in the shock comoving frame) are given by 



4TAGne 



(84) 



The electron distribution is assumed to be similar to that of internal shocks, N{'y) a 7 ^ for 
7m < 7 < 7Af, and we use p = 2.5 to obtain the numerical values. 

As mentioned in §4.4, the plerion radiation field is roughly homogeneous and isotropic within 

the the radius fRf, where the plerion emission takes place. At r < fR^ the external (plerion) photon 
energy density, Uph^ext, is given by equation (54), and we may use the relation U^f^^^^^ = F\G^ph,ext 
(that is valid for an isotropic radiation field) to obtain 



X 



TP 

ph,ext 

€Be' 



9 1 ^1/2 f- 2„-l -1/2 1/2 1 ^ a 



k = 



1 0^1/2 f-2 -4/3 -1/2 1/2 1 ^ p-1/3 ^2/3 -1/3 ^1/3 1/3 , _ -, 

k-L-O" /l/3''2/3 ^6e,l/3^6S,-3^S,-3Se,l/3-C'rot,53/^6,-l''s(i,1.5-^iso,53''days — J- 



1/3 ^2/3 ,-1/3 pl/3 ,1/3 



(85) 

where tsd,i.5 = isd/lO^'^- For tgd > ta ~ 65 yr the PWB is slow cooling, and / = 1, so that 
Rag < fRb throughout the afterglow. For tsd < ^3 the plerion emission is radiated within a shell 
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of width Rf, — Rg times the ratio, (isd/^s)^; of the coohng time of electrons with 75^, and the 
dynamical time t^d- This implies that generally, / = min{l, {Rs/Rb)[l — (^sd/^s)^] + (^sd/^s)^}- For 
tsd < ts we have / < 1 and therefore Rag < fRb only at sufficiently early times after the GRB. 
For tsd ^ ^3) the radiation is emitted within a thin shell behind the wind termination shock, at Rs, 
and / « Rg/Rb- In this case Rag > fRb throughout the afterglow. We study the implications in 
the following. 

At r > fRb, the plerion radiation field is no longer isotropic or homogeneous, and we model 
the plerion radiation field as resulting from emission by a uniformly bright sphere with a radius 
fRb, and obtain 



where = [1 — (/-Rft/r)^]^/^ and /3 = (1 — F^^)^/^. The ratio of the photon energy in the local and 
the observer frames is now 

( Fic(l - /^)V3 ~ ifRb/rfT\an2 r « Tag fRb 

■ (87) 

AG --^ -n - fT^ 

During the early afterglow, Rag is relatively small and Fag 3> 1, so that Rag ^ FAGfRhi and 
the first limit of equation (87) is applicable, implying 



%^ _ [(1 - ;3m)3 - (1 - /?)3] 

Uph 3/3(1 - /x) 



(1 - /?)/(! + 13)^ 1/4F2 g r » Tag fRb 



n Ofi„-l/2^4 1/2 -1/2 1/2 1 ^ p3/2 ^5/2 3/2 ^-3/2 -3/2 

U.^Iba /i/3'?2/3^6e,l/3^6S,-3^B,-3?e,l/3-f^rot,53P6,-l*sd,1.5^iso,53Ws ~ ^ 



nqo„-l/2M 2/3-1/2 1/2 1 ^ p5/3 ^8/3 5/3 ^-5/3 ,-5/3 
U.y<^a /l/3'?2/3Sea/3^6i?-3^B -3?e,l/3^rot,53/^b-lVl.5-^iso,53*days " ^ 

During the course of the afterglow its radius increases while its Lorentz factor decreases, so that 
eventually Rag becomes larger than TAGfRb-, and the second limit of equation (87) becomes 
relevant, implying 

' 1.3 X 10-^ a~^^^e^g^03ej^^_3e^^_34_l/3£^rot,53*3dji.5^iso,53idays ^ = 

^2 ~ < , (89) 

, 2.0 X 10-^ a~^'^^e^e!^/3^&B!-3eB^3^e,l/3-E^rot,53^Jd^l.5£^i;o,53^days ^ = 1 

as long as the afterglow shock is still relativistic. One can combine the two limits and use X = 
max(Xi,X2). However, since the region where these asymptotic expressions for X are not a very 
good approximation may play an important role, we use equation (86) rather than equations (88) 
and (89) for all our calculations. 

It is also worth to note that the average shift in frequency of the photons between the observer 
frame and local and rest frame is (y' jv) = FyiG'[l — /9(l + ;u)/2] [1 + T\q{1 — fi)'^]/2rAG, and varies 
between FyiG/2 and l/2FyiG- This should be compared to the usual factor of Tag for an isotropic 
(plerion) radiation field, and implies lower typical EC frequencies, by a factor of [1 — /3{1 + 1^)/2]. 
For simplicity we do not include this factor in the expressions for the EC frequencies, but we do 
take it into account in Figure 2, and when deriving constraints on the model parameters. 
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The electron cooling time is tsyn/ {l+Y+X) where the Compton y-parameter may be obtained 
by solving the equation 



aep 



eeil + Y + X) ' 



(90) 



which gives (Granot &; Konigl 2001): 



Y 



i/aee/es l,X^<aee/eB 
aeg/eB X, aee/f-B ^ 1 
^ oee/(esX) aee/es, K 



(91) 



For our choice of parameters, es <C eg and aee/es ^ 1, so that either the first or the third limits 
of Eq. (91) are relevant. As the first limit is more often applicable, we use the parameterization 
(1 + y + X) = A'(a€e/es)^/^, so that the numerical coefficients and explicit dependence on the 
parameters of the break frequencies (that depend on the electron cooling time), would be relevant 
for X'^ <^ aee/eB, where X 1. In the limit X'^ S> aee/eB, the numerical coefficients and parameter 
dependences change because of the dependence on X ^ X/(aee/eB)^/^, which is no longer close to 
1 in this limit. 



6.1. The Synchrotron Emission 

In the standard case of a uniform ambient medium, one can express the break frequencies and 
the peak flux in terms of the shock energy E, the ambient density nn, the observed time t, as well as 
€e, cb, and the distance to the source (Sari, Piran & Narayan 1998). This is thanks to the fact that 
for a 'standard' external medium that is composed of equal numbers of protons and electrons, so 
that both the shock dynamics (that is determined by li; = njiiJipCp) and the external number density 
of electrons (that enter the expressions for the flux normalization and self absorption frequency), 
are determined by a single parameter, nn- In the case of a shock propagating inside a PWB, the 
dynamics of the AG shock are determined hy w = ?T.H,cquiv"T^pC^ ~ Vlutn-pmpC^ that is dominated 
by the internal energy of the hot protons, while the number density of electrons is different, and 
dominated by the electron-positron pairs. We find 

Tj ^Yot, 3 3 

^^e'lTT'p ^ p. -1 ^ -1 rf^r,\ 

;^;r~ " I^TT^ '^2/3ee,l/37^,4.5 • (92) 

"■Hjequiv ^'/e Jw'"'e 

For an elongated PWB we can make the usual substitution = 3i?eff/2ctsd) to obtain the relevant 
expressions (see discussion just before §4.1). 

The self absorption frequency is typically Vsa < niin(^'c, z^m); and is calculated using equation 
(30) for fast cooling and equation (34) for slow cooling, solving for t[, = a' Rag AG = 1 for i^ia 
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and then Uga = ^AG^'sa- transition time from fast to slow cooling, Iq, is obtained by equating 
7c and 7^. For A; = we get 

V, ^ 2.1 X 10^^ (1 + z)-'/\-'X-\-l^e:^_,e-^]LlE;^^^^^^^^ Hz , 

^ 1.2 X 10^^ (1 + ^)'/%2'/3ee:i%e^,-l4(-37' ,4.54i'53C^ ' 

~ 9 V in21 n -U r^-5/8„-l/2i^-l„-V8 -1/2 1/2 p-1/8 ^3/8 3/8 ^1/8 ^-3/8 „ 
VM ~ 5.2 X lU (\ + Z) 'a ' X r?2/3 e^_^ -3^rot,53/^b -l*sd,1.5^iso,53*days ' 

. 1.0xlO«(l + .)-V2X,VjeS^e/^l^K"<"^^^ 

..„2 - 2.6xl0^(l + z)-S,73^;/f/3e-i,e^/f_3<,37-l'^^^^^^^^^ (93) 
i^,.,me« ~ 0.61 (1 + z)r]~ll'^^^^y^e^'^_^E]ll^^^-\^^^^ , 

to ~ 5.7 (1 + ^)X27y3^3e-2^3e^__,eB,_3^rot,537^,4.5^67!l*s11l.5^iso,53 days , 

where z/gai is for fast cooling and ^502 is for slow cooling. For k=l we get 

,/ ~ 70 V inl3 n -U r^-5/6„-l:^^-2 -4/3 1 -1/2 p-4/3 ^8/3 8/3 -1/6 -1/6 

Uc ~ 7.U X lU [I + Z) 'a X 772/3 «e -l^B -3^rot,53/^6 -l*sd,1.5'^iso,53*days ' 

« 9.9 X 10^4 (1 + ^)^/'^2'/3^e":i/3^e,-l4(-37^,4.54&3Cs' ' 

~ 4 9 V in21 n 4- r^-2/3„-l/2i^-l -1/6 -1/2 1/2 ^,-1/6 ^1/3 1/3 ^1/6 .-1/3 tt^ 

,y , ~ X inSri +rr2A5v„WlV3/5 ,V2 7/10 p22/15 -3/5^-44/15 -44/15 p^^^ -13/15 „ 

Z^sal ~ d.d X iU (i + 2;j -^^2/3 ^e,l/3^e -l^B -3-^rot,537«,,4.5P6 -1 *sd,1.5 -^iso,53*days ' 

, « 5 2 X 10^ ri + zY^l^n'^l^^^l^ e^^^ ^^/^ .,-8/5^-8/5 -8/5 -1/5 ^ 

Vsa2 ~ O.Z X iU \^^Z) r\^i^ '?e,l/3^e -l^B ,-3-^rot,537«,,4.5Pfe -1 '^sd,1.5''days ' l^^/* 

^ 1.2 (1 + z)'lS-,]i\.,i/^^l%E]ll,,^^^^^^^^ mJy , 

to - 7.3 (1 + .)X3/^4^e;J/Je^/^,4/;_3£;,„,,37:a5/?.!l^s-d!l.5<53 days , 

We note that in order for v^a not to exceed a few GHz, as typically implied by observations, we 
need tgd ^ 10 yr. This also gives more reasonable values for the transition time from fast to slow 
cooling, and for Fj^^max- For /c = the effective mass density of the PWB, becomes similar to 
that of a typical ISM, nH,equiv = 1 cm^^, for 

tiSM = 38 7?2/3^roi3/56;-l , (95) 

while the electron number density reaches the same value for a smaller tg^ = ^^•'^ yr- 
For tisM the afterglow emission is close to that of the 'standard' model, where the external 

medium is the ISM or a stellar wind, which has been extensively and successfully fitted to afterglow 
observations. 

In order to explain the X-ray lines, we need t^^ ^ ^Fc ~ 1 yr- This implies that the radio will 
typically be below the self absorption frequency, and hence the radio emission from the afterglow 
would not be detectable. On top of this, the jet break time is given by substituting nH,equiv in place 
of the external density for a 'standard' external medium (Sari, Piran h Halpern 1999), 

V 2 / \nH,equiv,0/ 
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0.70 



1 + z 



iso,53 



1/3 



2 J V%/3-^rot,53, 

and is very low for 9j-i = 6j/0.1 ~ 1. If we want to explain the observed values of tj ~ 1 day that 
are typically observed as resulting from a larger 9j (~ 0.4) then this would imply that the time of 
transition to a non-relativistic flow should be iNR ~ ^J^^j ~ 7tj ~ 7 days, and in general, 



E 



(97) 



^ \ '^H^equiv'^p^ 

where E = 10^^ £^51 ergs is the true energy of the afterglow, and we have dropped factors of order 
unity. Finally, for tgd ipe ~ 1 yr, the transition time from fast to slow cooling is very large, and 
fast cooling is expected during all the afterglow. 

For t <to, in the fast cooling regime, the synchrotron flux density, F^,, is given by-*^^: 



Ell — -^'i/,max ^ * 



(z./z.e)V3 



(98) 



For t > to we are in the slow cooling regime, in this case the spectrum peaks at and again 
consists of four power law segments: 



Ev — -^fc'.max ^ 



(99) 



^^If there is no significant mixing of the shocked fluid the the spectral slope just below v^a should be v^^l^ , and 
the familiar slope is obtained below a lower break frequency, Vac (4). 
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6.2. The SSC Emission 



The SSC emission is calculated similarly to §4.2 and §4.3. The fast cooling spectrum is given 



by 



•sc 



Y X I 



sa J 



SC 



u < v. 



(100) 



where 



VrnF.^ = 3.1 X 10-^^ (1 + ^)X-'eef-ies,-3^iso,53idays^Z2'8erg Cm"^ S 

and for = we have 

3.9 X 10^^ 



-2 „-l 



(101) 



.sc 



.sc 



.sc 



:^-lr,-3/4t3/5 -1/2 -3/10 p-3/20 -3/5^9/20 9/20 p-1/20 -1/4 „ 

~ (1 + ^)3/4 ''2/3 ?e,l/3^e,-l -3 -^rot,53 7«,,4.5A^6,-l'^sd,l-5-^iso,53 Mays 

^ IQ^^ „-2 ^-4„-9/4 2 ^-3/2 ^-9/4 ^27/4 27/4 -5/4 -1/4 „ 
(l + z)3/^ ^e,-l^B -3-^rot,53P6,-l'^sd,1.5'^iso,53Mays ' 



^,2„ 5.1 X 10^^ ^15/4^,4 4 1/2 p-1/4 4 ^3/4 3/4 3/4 ,-9/4 , . 

(I+^j-V^ ^2/3 ?e,l/3^e-ies -3-^rot,537«;,4.5P6-lIsd,1.5'^iso,53'^days -"^^ '."'■^^J 



sc 

^KN,1 



^AG^KNfi _ 6.7X1026 _9/4 



-z -l/2p-l/4 2 /,3/4 3/4 -x/^o/^ „ 
?e,l/3^e,-l^B,-3^rot,537«,,4.5/^fe,-l''sd,1.5-^iso,53Mays ' 



1/4,3/4 



Um (1 + 2) 



while for = 1 we have 



sc 



.sc 



^KNA 



3.3 X 10^^ , _4/5^3/5 1/2 -3/10^-1/5 -3/5^2/5 2/5 .-1/5 „ 

(1 _|_ 2;)4/5 '2/3 ^e,l/3^e,-l ^5,-3 -^rot,537«,,4.5P6,-l^d,1.5''days ' 

1.3 X 10^1 2 -_4 3 2 -3/2 p-3 ^6 .6 p-1/2,1/2 „ 

(1 + z)3/2 ^^2/3^6 -l^S -3-^rot,53/^6 -l^sd.l.S-f^iso.SS^days ' 

2.7 X 10^3 11/3^-4 4 1/2 p-1/3 4 ^2/3 2/3 ^5/6 .-13/6 „ . ^ 

(1 + z)-"^/^ 2/3 'Se,l/3^e -l^B -3-^rot,537«),4.5P6 -l'^sd,1.5-^iso,53'^days vJ-U-a; 

5.3 X 1026 -7/3^2 -2 -1/2 -1/3 -2 ^2/3 2/3 p-1/6 5/6 „ 

(1 + 2:)^^/6 2/3 ?e,l/3^e,-l%,-3-^rot,537i(;,4.5M6,-lSd,1.5-^iso,53''days ' 

2.5 X 1028 ^_i/4- ,1/2^-5/12,-1/4^^-5/12 ^5/6 ^5/6 pl/6 -1/12 „ 

(1 + Z)1V12 '2/3 e,-l -^rot,53 ,-l''sd,L5-'^iso,53''days ""^ " 
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The slow cooling spectrum is 

(,SC/,5C)(3-p)/2(,/,5C)4/3 
(z//z.f^)(3-P)/2 



ypSC 



sc 



V < V 



sc <; <- i^SC 



SC ^ <- jySC 



(104) 



where fc-^l/c is just a = (i^m/'^c)*-^"^^^^ times I'm.Fum for the fast cooling, that is given in Eq. (101). 
For A; = we have 



.SC 



sc 

^KNA 



1.1 X IQi^ 



_ ^3/4^-2/5, 1/5 p7/20 2/5 ^-21/20 -21/20 p9/20 -3/4 „ 

7m^sa2 ~ ^-^ _|_ ^^1/4 '/2/3?e,l/3^e,-ltB _3-C/j.ot,537u,,4.5P6,-l '^sd.l.S -^180,53^378 

riG4iV,0 3.8 X 1027 



.,^2 3/4 1/2 3/4 ^-9/4 -9/4 3/4 .-1/4 „ (.n.r■^ 

(1+^)3/4 ^2/3^e -l^S -3-^rot,53P6 -1 *sd,1.5 A80,53*days ' ^^^^^ 



and for A; = 1 we obtain 



.SC 



L.4 X 10^^ 13/15^-2/5 1/5 p7/15 2/5 ^-14/15 ,-14/15 pl/3 .-13/15 „ 

1 _). 2;)2/15 ^2/3 ^e,l/3^e'-l^B,-3-^rot,537«,,4.5^'6,-l ''8d,1.5 -^i80,53''day8 



( 

7.6 X 1027 



^,SC _ '-"^^^ „y2^ ^ ^1/2 p ^-2 ,-2 pl/2 -1/^ 
^i^Af.l ~ (1 +z)V2 ^2/3^6 -ies_3£'rot,53Pb_i«sd,1.5-^i80,53*day8 



3-2 ^-2 f1/2 ,-1/2 



Hz 



(106) 



6.3. The EC Emission 

The EC emission in this case arises from the upscattcring of the plcrion radiation by the 
relativistic electrons behind the afterglow shock. We provide detailed expressions for one repre- 
sentative plerion spectrum, the one given in equation (24). This spectrum is the spectrum for 
12 yr ~ t2 < i8d < ~ 65 yr (see equations 35, 23), which is of most interest. Similar expressions 
for the other plerion spectra can be readily derived in a similar manner. We note that for t^d < h 
the synchrotron emission of the plerion near the peak of vF^ is the same as for the spectrum we 
use (i.e. for t2 < tsd < ^3); a-iid therefore the EC near the peak of its I'Fi, should be the same. The 
peak of vFi, for the SSC plerion emission is typically above the KN limit for the AG electrons, and 
should therefore have a negligible contribution for the EC emission of the afterglow. The resulting 
EC spectrum depends on whether the electrons in the afterglow shock are in the fast cooling or 
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slow cooling regime. The EC spectrum is 



vF, 



EC 



EC 



= X X I 



u < v. 



(107) 



where ^'^^(tm) = ^^xatItm) is given by equation (102). For fast coohng and k=0 

iQ tV3 1/6 o23/6,7/2 

2.4 X 10^^ '^e,l/3^6B -3/^6,-1 *sd,1.5 



EC 



.EC 



EC 



^'^S4l/3^e.-l^S,-3-Ej/^537^(L4(^53*iays 



Hz 



6.4 X 10^9 



^'fe,-l''sd,1.5 



V^in X2^2%ef>e,l/3eb£_3ee,-ieiJ,-3^rlt,53^io%*yay; 



Hz , (108) 



6.5 X 10^1 



i AG7m^&m ~ ■^^-^P^jZT7^%/3Se,l/3^6e,l/3^6B,-3^e-l7ty,4.5-^^iso,53Ws 

while for fast cooling with A: = 1 we have 
..EC , 



1/2 



days 
,1/2 ,-3/2 



4.1 X 10^^ 2^1/3 -1/6 1/6 ,-1.-1 p;-5/3 -1/3^10/3 3 tTz 

^ ''2/3^e,l/3^be,l/3^ftS,-3^e,-l^S,-3-^rot,537io, 4.5/^6,-1 '-sda-S -"^ ' 



EC 



.EC 



(1 + ^) 

1.1 X 10^9 



1/2 -1 



e-l e-l ^-7/2 ^17/2 13/2 
^e,-l^B,-3-^rot,53'^6,-l''sd,l-5 ' 



<r-2 -1 

+ 2) ^2/3 ^6e,l/3^6B,-3 

2.7X10^1 23/6^-4 .2 ^1/2 2 ^"1/6 4 ^-5/6 -5/6^2/3 .-4/3 „ 

,l/3^6B,-3^e,-l-^rot,53 '«),4.5^'6,-l ''sd,l-5-^iso,53''days 



(109) 



(l + z)-V3 ^2/3 ?e,l/3^6e 

For slow cooling with A; = we have 



.EC 



.EC 



.EC 



2 2 
^AGlc^bm 



-.r.^Q J/^.^/'^ .2 5/3 ^5/6 1/2 1/2 
, 1.4 X 10^9 ^2/3^bB -3ge,-l7^»,4.5Pb,-l^sd,1.5-^iso,53 

' fl + z)-V2 ;^5/3 1/6 pl/6 3/2 

^ ^ ?e,l/3^6e,l/3-^rot,53''days 

3.7x1020 fr,-i;ii.5t-i'ii,i.5^i!^53 



Hz 



Hz , 



(1 + Z)-V2 /:2 , , ,1/2 p2 ^3/2 
^ ' ?e,l/3^6e,l/3^6B-3'^rot,53^days 

1.1 X 10^1 ^2/3eL,l/3eb£-37^,4.5/5fe,-l*sd,1.5 



(110) 



+ ^ aX2^2^^_^€e -les -3^rot,53^ioi3*iays 



Hz , 



and finally for A; = 1 we have 



.EC 



X 10^^ 4/3^-5/3 -1/6 1/6 2 p-1/3 5/3 ^2/3 1/3 ^2/3 -4/3 
(1 + z)-^/^ ''2/3?e,l/3^6e,l/3^6B,-3^e,-l-^rot,537u,,4.5A^6,-l''sd,1.5-^iso,53''days 
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,EC 



EC 



-1/6^-2 ^-1 ^-1/2 2 p-13/6 2 ^35/6 23/6 ^2/3 .-4/3 ^r^ 
''2/3 ^e,l/3^6e,l/3^6B,-3^e,-l-^rot,53 7m),4.5P6 -l'^sd,1.5-^iso,53''days ^V^^^) 



-1/2 2 



,-13/6 2 



1.6 X lO^o 

(1 + z)-V3 ^2/3 ^e,l/3'=6e,l/3'=6B,-3^e,-l--rot 

1.9 X 10^° _i - _2 1/2^-2 2 1/2 -1 -1 p-3/2 2 ^35/2 5/2 „ 

^ ^ '^2/3^e,l/3^&e,l/3^6B,-3^e,-l^S,-3-^rot,537«),4.5P6,-l'^sd,1.5 



(1 + ^) 



6.4. High Energy Emission 

Figure 2 shows the uFi, spectrum of the afterglow at t = 500 s and 5 x 10^ s, for our fiducial 
parameters, and for tgd = 20 yr, £'rot,53 = 0.5, z = 1, Rg/Rb = 0.1. As can be seen from the figure, 
for t = 500 s (5 X 10^ s) the synchrotron is dominant below 45 MeV (2 MeV), while the EC is 
dominant above this range. At t = 5 x 10^ s the SSC component becomes dominant above 40 GcV. 
We expect an upper cutoff at hv^y^ ~ 250 GeV, due to opacity to pair production, with the photons 
of the plerion. This upper cutoff moves down to a lower energy for smaller values of tsd, and is 
~ 1 GeV for tsd = 1 yr ~ ipe- This implies that for afterglows with X-ray line features we expect 
no high energy emission above this limit. 

We find that the early afterglow (t < 100 s) emission at > 100 MeV is dominated by the 
EC and SSC component, which are comparable at this time. At later times the EC becomes 

dominant over the SSC component. The peak of the vF^'-' emission is at the level of ~ 5 x 
10"^(t/500 s)-^ ergscm-2s"\ and is located at hv^^ ~ 70(^/500 s)-3/2 GeV (see Eq. 108). 
The spectrum scales as vF^ oc jy^^"*)/^ above this photon energy, implying a flat vFy for values 
of s 2 that are typically inferred for PWBs, while it scales as vFy oc v^^"^ below this energy. 
At early times the afterglow radius is relatively small (-Rag ^ fRb)i so that X is approximately 
constant in time, and the peak of the vF^, EC spectrum has a temporal scaling similar to that for 
the synchrotron component (i.e. oc see Eq. 101). We therefore expect vFjj at a fixed photon 
energy to decay very slowly with time, as at v < ly^f" , and decay approximately linearly with 

time (oc i~^~3(s-2)/4-j > t'm*^- The temporal decay becomes larger than these scalings as the 
afterglow radius Rag increases above fRb, and the parameter X begins to decrease with time. 

This result can explain the high energy emission detected by EGRET for GRB 940217 (Hurley 
et al. 1994). This detection consists of photons of energies > 50 MeV, where 10 photons were 
observed during the prompt GRB emission, that lasted 180 s, and 18 photons were detected up 
to 5, 400 s after the end of the GRB, which include a photon of energy 18 GeV. During the post- 
GRB emission, the source position was Earth-occulted for ~ 3, 700 s. At t ~ 500 s the flux is 
~ 1 — 2 X 10~^ ergs cm~^ s~^ which is roughly consistent with our results. At t 5, 000 s, after the 
Earth-occultation, the flux is a factor of ~ 2 — 3 lower, if we exclude the one 18 GeV photon. This 
moderate time decay is consistent with our results. 

A different interpretation for the high energy emission discussed above was recently suggested 
by Wang, Dai &: Lu (2002), in a similar context of the supranova model, where the GRB occurs 
inside a plerionic environment. The main difference is that they consider a pulsar wind that consists 
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purely of e pairs, that is in the fast cooling regime, and therefore the radius of the termination 
shock of the wind, Rg is very close to the outer radius of the PWB, Rh, and all the shocked wind 
is concentrated within a thin radial interval, Rh — Rg <^ Rf,. They try to explain the high energy 
emission as the synchrotron emission from the early afterglow. They obtain hi'rn ~ 1 GeV at 
t ~ 180 s, and according to their model t ~ Ri,/2T'^c or F oc t"^/'^ and oc t""^. However, 

this implies hv^ ~ a few keV after one day, which is inconsistent with afterglow observations. They 
also claim that the EC emission is unimportant, which is in contradiction with our conclusions. 
The inclusion of a proton component in the pulsar wind with a similar energy to that of the 
component allows only the energy in the to be radiated away, so that even for a fast cooling 
PWB a large part of the energy of the pulsar wind remains in the protons and Rg is significantly 
smaller than R^. 



7. Discussion 

We have studied the observational implications of GRBs occurring inside a pulsar wind bubble 
(PWB), as is expected in the supranova model. We find that the most important parameter that 

determines the behavior of the system is the time delay, tgdj between the supernova and GRB events. 
The value of tgd is given by the typical timcscalc on which the SMNS loses its rotational energy 
due to magnetic dipolc radiation (sec Eq. 2) and depends mainly on the magnetic field strength of 
the SMNS, (since its mass, radius and spin period are constrained to a much smaller range of 
possible values). For ~ 10^^ — 10^^ G, t^A is between a few weeks and several years. However, a 
larger range in B*, and correspondingly in tgd) seems plausible. We therefore consider t^d as a free 
parameter, and predict the observational consequences of different values for this parameter: 

1. For extremely small values of tgd < ^coi = R^/PbC ~ 0.9 Ri,^i3/3^\ hr, where i?* = 
lO^^R^^is cm is the radius of the progenitor star (before it explodes in a supernova), the stel- 
lar envelope does not have enough time to increase its radius considerably before the GRB goes 
off, and the supranova model reduces to the coUapsar model. In this respect, the collapsar model 
may be seen as a special case of the supranova model. Such low values of tgd niight be achieved if 
the SMNS is not rotating uniformly, as differential rotation may amplify the magnetic field to very 
large values. 

2. When tcoi < ^sd < ^IS ~ 3/3^^i^r^^ « l6P^^\rl^^ty-2 days (for k = 0) the deceleration 
radius i?dec ~ i^NRP"^''^^"'^^ ~ i^bF^^/^^"*^^ is smaller than the radius for internal shocks Rjs ~ 
2F^ct„. In this case the kinetic energy of the GRB ejecta is dissipated through an external shock 
that is driven into the shocked pulsar wind, before internal shocks that result from variability within 

the outflow have time to occur. For an elongated PWB, tjs can be smaller by up to a factor for 
~ 10, since the polar radius would be ~ 10 time larger for the same tgd, and the volume of the 
PWB would be much larger, and the density much smaller. 

3. If tis < tsd < ~ 0.4 yr, internal shock can occur inside the PWB, but the SNR shell is still 
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optically thick to Thomson scattering, and the radiation from the plerion, the prompt GRB and 
the afterglow cannot escape and reach the observer. If the SNR shell is clumpy (possibly due to 
the Rayleigh- Taylor instability, see §2), then the Thomson optical depth in the under-dense regions 
within the shell may decrease below unity at tgd somewhat smaller than tr, enabling some of the 
radiation from the plerion to escape. For an elongated PWB, the polar radius can be larger by up 
to a factor of ~ 10, which reduces tr by the same factor. Furthermore, the elongation can be due 
to a smaller than average surface mass density of the SNR shell at the poles. This would further 
reduce tr- 

4. For tr < tgd < tpe ~ 1 yr the SNR shell has a Thomson optical depth smaller than unity, 
but the optical depth for the iron line features is still > 1 so that detectable X-ray line features, 
like the iron lines observed in several afterglows, can be produced. 

5. Finally, for tpc < ^sd) we expect no iron lines. When tg^ is between ~ 2 yr and ~ 20 yr the 
radio emission of the plerion may be detectable for 7^ < 10^. The lack of detection of such a radio 
emission excludes values of igd in t^iis range, if indeed 7^, < 10^, as is needed to obtain reasonable 
values for the break frequencies of the afterglow. For tgd = iiSM ~ 38 yr, the effective density of 
the PWB is similar to that of the ISM (i.e. 1 cm~^), and the afterglow emission is similar to that of 
the standard model, where A; = is similar to an ISM environment, with the exception that in our 
model a value of A; = 1, that is intermediate between an ISM and a stellar wind, is also possible. 
Larger (smaller) values of the external density are obtained for smaller (larger) values of tgd- 

The SNR shell is decelerated due to the sweeping up of the surrounding medium for tgd > 
225 Mg^'^ ^^"^/5^1]^ yr, where n = uq cm^^ is the number density of the external medium, which 
is larger than the values of igd that are of interest to us. This effect may therefore be neglected for 
our purposes. 

An important difference between our analysis and previous works (KG; IGP) is that we allow for 
a proton component in the pulsar wind, that carries a significant fraction of its energy. In contrast 
to the 6=*= component, the internal energy of the protons in the shocked wind is not radiated away, 
and therefore a large fraction of the energy in the pulsar wind (~ 10^^ ergs) is always left in the 
PWB. This implies that even for a fast cooling PWB, the radius of the wind termination shock 
Rs is significantly smaller than the radius of the SNR, Rf,, and that the afterglow shock typically 
becomes non-relativistic before it reaches the outer boundary of the PWB. 

8. Conclusions 

Our main conclusion is that existing afterglow observations put interesting constraints on the 
model parameters, the most important of which being the time delay igd between the supernova 
and GRB events, which is constrained to be > 20 yr, in order to explain typical afterglow observa- 
tions and the lack of detection of the plerion emission in the radio during the afterglow. Another 
important conclusion is that iron line features that have been observed in a few X-ray afterglows 
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cannot be naturally explained within the simplest spherical version of the PWB model, that has 
been considered in this work. This is because the production of these lines requires ^ ^Fc ~ 1 yr 
which implies a very large density for the PWB and effects the afterglow emission in a number of 
different ways: i) The self absorption frequency of the afterglow is typically above the radio, imply- 
ing no detectable radio afterglow, while radio afterglows were detected for GRBs 970508, 970828, 
and 991216, where the iron line feature for the latest of these three is the most significant detection 
to date (~ 4(7). We also expect the self absorption frequency of the plerion emission to be above 
the radio in this case, so that the radio emission from the plerion should not be detectable, and 
possibly confused with that of the afterglow, ii) A short jet break time tj and a relatively short 
non-relativistic transition time tjsiR are implied, as both scale linearly with igd and are in the right 
range inferred from observations for tgd ~ 30 yr (see Eqs. 96, 97). iii) The electrons are always in 
the fast cooling regime during the entire afterglow. 

The above constraints regarding the iron lines may be relaxed if we allow for deviations from 
the simple spherical geometry we have assumed for the PWB. A natural variant is when the PWB 
becomes elongated along its rotational axis (KG). This may occur if the surface mass density of 
the SNR shell is smaller at the poles compared to the equator, so that during the acceleration of 
the SNR shell by the pressure of the shocked pulsar wind (that is expected to be roughly the same 
at the poles and at the equator) its radius will become larger at the poles, as the acceleration there 
will be larger. A large-scale toroidal magnetic field within the PWB may also contribute to the 
elongation of the SNR shell along its polar axis (KG). It is also likely that the progenitor star that 
gave rise to a SMNS had an anisotropic mass loss, which results in a density contrast between the 
equators (where the density is higher) and the poles (where the density is lower). A sufficiently 
large density contrast between the equator and the poles can also contribute to the elongation of the 
shell, for sufficiently large tgd, as the SNR shell will begin to be decelerated due to the interaction 
with the external medium, at a smaller radius near the equator, compared to the poles. A similar 
non-spherical variant of the model is if we allow for holes in the SNR shell, that extend over a small 
angle around the polar axis, where all the wind is decelerated in a termination shock within the 
SNR shell (Rg < Rb), but most of the shocked pulsar wind can get out through the holes near the 
poles and reach a radius considerably larger than i?^. This variant may be viewed as a limiting case 
of the previous variant, when the surface density contrast of the SNR shell, between the equator 
and the poles is very large. This implies that most of the mass in the SNR shell is concentrated 
near the equator, while only a small fraction of it is near the poles, so that the radius near the poles 
can be as large as ~ ctsd, while the equatorial radius is < R^. In both variants, the total volume of 
the PWB is close to that of a sphere with the polar radius, and much larger than that of a sphere 
with the equatorial radius. This would allow for a smaller density with the same small value of 
the equatorial radius that is required to produce the iron lines. We see that in principle, variants 
of the simple model are capable of reconciling between the iron line detections and the afterglow 
observations. 

An important advantage of the PWB model is that it can naturally explain the large values of 
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Cfi and eg that are inferred from fits to afterglow data (KG), thanks to the large relative number 
of electron-positron pairs and large magnetic fields in the PWB. This is in contrast with standard 
environment that is usually assumed to be either an ISM or the stellar wind of a massive progenitor, 
that consists of protons and electrons, and where the magnetic field is too small to explain the values 
inferred from afterglow observations, and magnetic field generation at the shock itself is required. 
Additional advantages of the PWB model are its ability to naturally account for the range of 
external densities inferred from afterglow fits, and allowing for a homogeneous external medium 
{k ^ 0), as inferred for most afterglows, with a massive star progenitor. 

Another advantage of the PWB model is its capability of explaining the high energy emission 
observed in some GRBs (Schneid et al. 1992; Sommer et al. 1994; Hurley et al. 19994; Schneid 
et al. 1995). We find that the high energy emission during the early afterglow at photon energies 
> 100 keV is dominated by the external Compton (EC) component, that is due to the upscattering 
of photons from the plerion radiation field to higher energies by the relativistic electrons behind the 
afterglow shock. We predict that such a high energy emission may be detected in a large fraction 
of GRBs with the upcoming mission GLAST. However, we find an upper cutoff at a photon energy 
of 1 Q GeV, due to opacity to pair production with the photons of the PWB. This implies no 
high energy emission above ^ 1 GeV for afterglows with X-ray line features, but allows photons up 
to an energy of ~ 1 TeV for afterglows with an external density typical of the ISM (tg^ ~ iiSM)- 
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Fig. 1. — The flux density, Fj,, of the plerion emission at the time of the GRB (i.e. a time tgd after 
the supernova event), as a function of tgd- The three panels show the flux density in the radio, 
optical and X-ray bands (z/ = 5 x 10^, 5 x 10^^, 10^^ Hz, respectively). For the radio and optical 
we show Fy of the synchrotron emission, and for the X-ray we show vFy for both the synchrotron 
and SSC components. The dashed, solid and dot-dashed lines represent log^Q(7,„) = 4, 4.5 and 5, 
respectively. Dashed lines are shown at t^- (below which the Thomson optical depth is larger than 
unity), tFe (below which iron line features can appear in the X-ray spectrum of the afterglow) and 
iiSM (for which the effective density of the PWB is similar to that of a typical ISM, i.e. 1 cm~^. 



-43- 




10" 
10" 
810" 

(0 

10" 
-^10" 

D) _ 
felO 

.^10" 
LlIO" 

>10" 
10" 
10" 



"T — I — I 



ir| I I 



syn 
SSC 
EC 



t = 5 X 1 0^ sec 



■9 
■10 
■11 
■12 

■13 
■14 
■15 
■16 

10^^10^^0^^0^^10^^0^^0^^0^^0^°10^^10^^10^^0^^0^^10^^ 



yy 



ill I I 



V in Hz 



Fig. 2. — The afterglow uF^, spectrum at t = 500 s (upper panel) and 5 x 10^ s (lower panel) after 
the GRB, for our fiducial parameters, and for tgd = 20 yr, -Erot,53 = 0.5, z = 1, Rg/Rh = 0.1. 
Dotted vertical lines indicate um where the upper cutoff for the synchrotron emission is located, 
and Vry^ where the upper cutoff of the SSC and EC (due to pair opacity) is located. 



